Bethe Eigenvectors of Higher Transfer Matrices 



E.Mukhin*, V. Tarasov*, and A.Varchenko 



*>* Department of Mathematical Sciences 
Indiana University Purdue University at Indianapolis 
Indianapolis, IN 46202, USA 

*St. Petersburg Branch of Steklov Mathematical Institute 
Fontanka 27, St. Petersburg 191023, Russia 

^Department of Mathematics, University of North Carolina 
Chapel Hill, NC 27599, USA 

April 2006 

Abstract. We consider the XXX-type and Gaudin quantum integrable models as- 
sociated with the Lie algebra $l N . The models are defined on a tensor product 
Mi <g> . . . <g> M n of irreducible gl^-modules. For each model, there exist iV one- 
parameter families of commuting operators on Mi ® . . . ® M n , called the transfer 
matrices. We show that the Bethe vectors for these models, given by the algebraic 
nested Bethe ansatz, are eigenvectors of higher transfer matrices and compute the 
corresponding eigenvalues. 

1. Introduction 

In this paper we consider quantum integrable models associated with the Lie alge- 
bra $l N , the XXX-type and Gaudin models. The models are defined on a tensor prod- 
uct Mi ® . . . <g) M n of irreducible gt^-modules. For each model, there exist N one- 
parameter families of mutually commuting operators on Mi <g> . . . <g> M n , called the trans- 
fer matrices. In this paper the transfer matrices for the XXX-type model are denoted 
by 7^- min (u; z),..., 7^-® Mn (u; z) , see Section 6, and the transfer matrices for the 
Gaudin model are denoted by gMx®-®M n ( M; ^ ... ) gM±®...®M n ^. ^ ? gee g ectkm g 

The important problem is to find eigenvectors and eigenvalues of the transfer matrices. 
If Mi, . . . , M n are highest weight modules, this can be done by the algebraic Bethe ansatz 
method. The N = 2 case is fairly well understood, for example, see [G] , [KBI] . However, 
the N > 2 case is much more technically involved and less studied. 
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The purpose of this paper is to obtain eigenvectors and eigenvalues of the transfer ma- 
trices for N > 2 using the algebraic nested Bethe ansatz method. For the transfer matrix 

<jMi®...®M n £ u . z j thig wag done in j KR1 | _ transfer matrix jMi®...®M„ ^. ^ comeg from 

the central elements of the Yangian Y(gl N ) and is proportional to the identity operator on 
the whole tensor product Mi <g> . . . <g> M n [KS] . In this paper we extend the technique devel- 
oped in [KR1] to show that the Bethe vectors are eigenvectors of the higher transfer matrices 
fjfc^'"® Mn (u) , k > 1 , and to compute the eigenvalues, see Theorems 5.4 and 6.1. The for- 
mula for eigenvalues of the higher transfer matrices on Bethe vectors were conjectured in 
[KR2] . However, no proof can be found in literature. 

For the Gaudin model, the transfer matrix Q^j^'"® Mn (u] z) comes from the central el- 
ements of the current algebra gljvf^] an d is proportional to the identity operator on the 
whole tensor product Mi® . . . (g> M n . The transfer matrix Q^^'"® Mn (u; z) is the gener- 
ating function of the famous Gaudin Hamiltonians which also appear as the operators in the 
Knizhnik-Zamolodchikov equations. Eigenvectors and eigenvalues of Q^j^'"® Mn [u] z) can 

be obtained from those of 7^^'"® Mn {u; z) by taking a suitable limit. For more straight- 
forward approach see [Ju] . For a construction of eigenvectors and eigenvalues of the Gaudin 
Hamiltonians as a quasiclassical limit of solutions to the Knizhnik-Zamolodchikov equations 
see [RV] . 

The transfer matrices for the XXX-type model were constructed in [KS] more than 
20 years ago. Suprisingly, an explicit form of the higher transfer matrices for the Gaudin 
model became known only recently. In this paper we use the approach suggested in [T] 
that gives the transfer matrices of the Gaudin model as a suitable limit of proper linear 
combinations of the transfer matrices of the XXX-type model. We trace this limit in the 
proof of the Theorem 5.4. This allows us to show that the Bethe vectors of the Gaudin 
model are eigenvectors of the higher transfer matrices and to compute the eigenvalues; see 
Theorems 8.6 and 9.2. 

There is another approach to higher transfer matrices of the Gaudin model and the 
Bethe ansatz, based on the representation theory of affine Lie algebras at the critical level 
[FFR]. It is interesting and important to identify the transfer matrices Q^^'"® Mn {u- : z), 
. . . , gMi®—®M n z ^ amon g higher Gaudin Hamiltonians constructed in [FFR] . 

Following [T] , it is convenient to combine the transfer matrices of the XXX-type model 
into a single difference operator 

©J£®-® M " (u; d u ; z) = 1 - Tj£®-®^ ( u; z ) e ~ d - + ...+ (-1) N -® M " («; z) e~ Nd ™ , 

where d u = d/du, cf. (A.l), (4.16). The operator X>^§"^ Mn (u;d u ; z) acts on functions 
of u with values in M\ <£>...<£> M n . Then any common eigenvector v G M\ <8> . . . ® M n of 
the transfer matrices, 

0jg 8 -® M » (u;z)v =At(u)v, k = l,...,N, 

where A\{v), . . . , A v N {u) are suitable scalar functions, defines the scalar difference operator 

M v (u; d u ) = 1- A v i(u) e~ d - + ...+ (-l) N A v N (u) e~ Nd " . 
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If v = Mf®-® Vn (t;z) is the Bethe vector of the XXX-type model, cf. (6.2), (5.2), then 
the operator Wl v (u; d u ) coincides with the operator ^R^q{u, d u ; t; z; A) , cf. (6.6) , which is 
written as a product, 

(1.1) m^ Q (u,d u ;t;z;A) = (l - X\ tQ (u; t;z; A) e" a «) . . . (l - X^ Q («; t;z; A) e" 9 ") , 

where the functions X| q(w; t; z; A) , see (6.4) , are explicitly defined in terms of the solution 
t of the Bethe ansatz equations (6.3) . This gives a formula for the corresponding eigenvalues 
A\(u;z),...,A v N (u). 

It is also convenient to combine the transfer matrices of the Gaudin model into a single 
differential operator 

Df 8 - 8M " («; d u ; z) = - g™-®" n {u . z) Q*-i +mmm+ ( _ 1} iv g M^...®M n {u . z) ^ 
cf. (A. 2) , (8.3) . Then any common eigenvector v G Mi ® . . . <g) M n of the transfer matrices, 

= fc = l,...,X, 

where SJ'(u), . . . , 2^(u) are suitable scalar functions, defines the scalar differential operator 

m> ; ^) = af-^(n)af- 1 + ...+ (-i) Jv ^( w ). 

If w = Wf®-® Vn (t;z) is the Bethe vector of the Gaudin model, cf. (9.6), (8.7), then the 
operator M v (u;d u ) coincides with the operator M.£ : k(u, d u ; t; z; A) , cf. (9.4), which is 
written as a product, 

(1.2) Mt, K (u,d u ;t;z;A) = (d u - y \u; t; z; A)) . . . (d u - )) N (u; t; z; A)) , 

where the functions ]$ a (u;t; z; A) are explicitly defined in terms of the solution t of the 
Bethe ansatz equations (9.3) . This gives a formula for the eigenvalues Z\ (u), . . . , Z"^{u) . 

The Gaudin model can be obtained from the XXX-type model in the limit as «,z tend to 
infinity. In this limit, the difference operators D^ 1 Q'"^ lMn (u;d u ;z) and Wl^Q(u,d u ; t; z; A) 
turn into the differential operators £)^i®---® M ™ ^ anc [ jyj;^^ <9 U ; t; A) , 

S M 1 ®...®M„ (£ -l W;e ^. £ -l^ = £ N V Mi®...®M n ( u . du . z) + ( £ N+l^ 
m Mi®...®M n e _i z . A) = £ N M M^...®M„ (u; ^. t; z . A) + 0(^+1) j 

as Q = 1 + and £ — > . Moreover, the Bethe ansatz equations (6.3) of the XXX-type 
model become in this limit the Bethe ansatz equations (9.3) of the Gaudin model, and the 
Bethe vector B][ 1 ®-® t '"(t; z ) turns into the Bethe vector f^ 10 -® u " (£; z) , cf. Theorem A.l. 

The operators (1.1) and (1.2) of the XXX-type and Gaudin models appeared in [MV2] 
and [MV1] , respectively. Those operators were called in [MV1] , [MV2] the fundamental 



3 



operators of the corresponding solutions of the Bethe ansatz equations. It was shown therein 
that for finite-dimensional jjl^-modules Mi, . . . , M n , the kernel of a fundamental operator 
is generated by polynomials only and the singular points of the fundamental operator were 
described. These facts and Schubert calculus imply that the number of solutions of the 
corresponding Bethe ansatz equations (properly counted) does not exceed the dimension 
of Mi <g) ... <g) M n [MV1] , [MV2] . That was a step to the proof of the completeness of the 
Bethe ansatz conjecture, which says that the Bethe vectors form a basis of Mi <g> . . . ® M n 
under certain conditions. 

An important observation for applications is the fact that the transfer matrices are sym- 
metric with respect to a suitable symmetric bilinear form on Mi <g> . . . <g> M n , see Theo- 
rems 9.1 and 7.6. Under certain assumptions, this fact allows us to conclude that the 
eigenvalues of the transfer matrices take real values for real u . In particular, this implies 
that the polynomials generating the kernel of the fundamental operator of a solution of the 
Bethe ansatz equations can be chosen to have real coefficients only, see applications to real 
algebraic geometry in [MTV] . 

2. Basic notation 

We will be using the standard superscript notation for embeddings of tensor factors into 
tensor products. If Ai, . . . , Ak are unital associative algebras, and a G At , then 

a (0 = a 1 ®(fc-0 eAi®...®A k . 

If a e Ai and b G Aj , then (a <g> 6) = a® , etc. 

Example. Let k = 2 . Let Ai , A2 be two copies of the same algebra A . Then for any 
a,beA we have a (1) = a <g> 1 , 6 (2) = 1 <g) b, (a <g) 6) (12) = a <g> 6 and (a ® 6) (21) = 6 ® a. 

Fix a positive integer iV . All over the paper we use the convention V = C N , and we 
identify elements of End(U) with NxN matrices using the standard basis of C N . 
Let e a b , a, b = 1, . . . , N , be the standard generators of the Lie algebra gl N : 

[e a b , e c d] = 5 bc e ad - 5 ad e cb . 

A vector v in a gl^-module is a weight vector of weight (A 1 , . . . , A^) if e aa v = A a v for any 
a = 1, . . . , N. A vector v is called a singular vector if e a tf = for any 1 ^ a < b ^ N. 
Let A{ fc j. G EndfV* 8 ^) be the skew-symmetrization projector: 

A {k} v 1 . . . (8) ^ = — ^ sign(a) u <Tl <8> . . . <g> v CTfc . 

We denote by V Ak the image of A{ fc }, and for any Q e End(V) set Q Ak = Q m \ yAk - In 
particular, Q AN = det Q • 

The space V is considered as a gl A r" m °dule with the natural action: e a b ^ -^ab , where 
E a }, G End (V) is the matrix with the only nonzero entry equal to 1 at the intersection of 
the a-th row and 6-th column. The gl^-module V is called the vector representation. The 

space V Ak is a gl^-module with the action e a b 1— > (E^ + . . . + ^5^)1 Afc - 



4 



We will use the following notation for products of noncommuting factors: 



\ \ Xi = Xi . . . Xk , Y[ Xi = Xk . . . Xi , 

l^i^A; l^i^k 

n = n n n = n n 

l^*<j^fc l^i^k i<j^k l^iKj^k l^j^k l^i<j 



3. 12-matrices 



N 

Let P = Yl Eab <8> Eba • It is the flip map: P(x ®y)=y®x for any x, y G V, and 

a, 6=1 

^(21) = p A p for any ^ e End (y®2) _ 

The rational R-matrix is i?(/u) =m + Pg End(F® 2 ) . It satisfies the inversion relation 

(3.1) R(u)R {21 \-u) = 1-u 2 
and the Yang-Baxter equation 

(3.2) R^ 12 \u — v )R^ 13 \u)R^\v) = rW(v)R( 13 \u)R( 12 \ u — V . 

Remark. The .R-matrix R(u) is symmetric: R(u) = R( 21 \u) . However, we will not use 
this property and will write all relations, for instance (3.1), in the form that naturally 
extends to more general .R-matrices. 

Below we use the well known fusion procedure to define .R-matrices R Ah,Al (u) acting 
in the tensor products V Ak (g> V Al , and describe their properties. The proofs can be found 
in [MNO, Sections 1,2], [N, Sections 2,3]. Notice that our notation differs from those in 
[MNO] , [N] , instead of R(u) the .R-matrix 1 — P/u is used therein. 

Lemma 3.1. fj R^\i-j) = f[ R^\i - j) = A {k} (-l) fc J] (-j) k ~ j+1 . 

l^i<j^k l^i<j^k j=l 

The Yang-Baxter equation (3.2) and Lemma 3.1 imply that 

(3.3) Af k f ) Af+ 1 --- k+l) fl H R (iJ+k \u + i-j-k + l) = 

l^i^k l^j^l 

(l...k) (k+l...k+l) 



= ( n n r[ " ,+i) ^ +i-j-k+Q) A<i T x> + 

^ i<i<k ' 



LC7 

Set 

R Ak > Al (u) = fl R (t ' J+k Hu + t-j-k + l)\ vAk ^ vM eEnd(V Ak )®End(V Al ) 
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The introduced .R-matrices satisfy the inversion relation 

k I 

r^\u) (R^ k (-u)f i} = nn(i-(«-i+ j) 2 ) , 

i=l j=l 

and the Yang-Baxter equation 

(R Ak > Al (u - v)) {12) (R Ak > Am (u)) {13) (R Al > Am (v)) {23) = 

= (R Al > Am (v)) {23) (R Ak > Am (u)) {13) (R Ak > Al (u - v)) {12) . 
In addition, for any Q e End (V) , 

(3.4) [R Ak > Al (u),Q Ak ® Q Al ] = 0. 

Set 

JV 

J2a*,ai(«) = «+ E (^ ) +--- + ^2 ) )|^ fc ®^aGEnd(y Afc )®End(y), 

o,6=l 

and 

N 

RaiMu) =u + k-l+ £ K6®(^ ) +... + ^ ) )|^ Afc GEnd(y)®End(y Afc ). 

a, 6=1 

In particular, Ran,ai(u) = u + 1 and R ai, a n(u) = u + N . 

k-l k-2 

Lemma 3.2. # Afe ' Al (u) = i2 Afc , A i(u) J] (u - i) , # A1 ' Afc (w) = i2 A i, A fc(«) II ( w + z ) ■ 

i=l i=0 

Corollary 3.3. # A fc,Ai(u) (-R A i, A /c(-«)) = (it + 1) (fc - it) ■ 

4. Yangian FQjtjv) 

In this section we collect some known facts from the representation theory of the Yangian 
Y(qI n ) that will be used in the paper. We refer the reader to reviews [MNO] and [Mo] for 
proofs and details. Notice that the series T ab (u) in (4.1) corresponds to the series T ba (u) 
in [MNO] , [Mo] . 

The Yangian Y(qI n ) is the unital associative algebra with generators , a, b = 1, ... , 
N and s = 1, 2, . . . . Organize them into generating series: 

oo 

(4.1) T ah (u) =^+^T { 6 V, a, 6=1,..., N. 

s=l 

The defining relations in Y(gl N )) have the form 

(4.2) {u-v)[T ab (u),T cd {v)] = T cb (v)T ad (u)-T cb (u)T ad (v), 
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for all a,b,c,d = 1, . . . , N . 

N 

Combine all series (4.1) together into a series T{u) = Yl Eab <8> T a b(u) with coefficients 

a, 6=1 

in End(C Ar ) <E> Y(gl N ) . Relations (4.2) can be written as the following equality for series 
with coefficients in End(C 7V ) <g> End(C N ) <g> Y(gl N ) : 

(4.3) R {12) [u — v (u — V) . 

The Yangian Y(gl N ) is a Hopf algebra. In terms of generating series (4.1), the coproduct 
A : Y(gl N ) — »■ Y(gl N ) <S> Y(gl N ) reads as follows: 

N 

(4.4) A(T ab (u)) = T cb (u)®T ac (u), a, b = 1, . . . , N . 

c=l 

This formula amounts to an equality 

(4.5) (id <g) A) (T(u)) = T( 13 ) (u) T( 12 ) (u) 

for series with coefficients in End(F) ® ^(gtAr) <8 ^(flfjv) • 

There is a one-parameter family of automorphisms p x : Y(gl N ) — > y(flfjv) that is defined 
in terms of the series T(u) by the rule p x (T(u)) = T(u — x); in the right hand side, 
(u — x)~ l has to be expanded as a power series in -u -1 . 

The Yangian Y(gl N ) contains the universal enveloping algebra U(gl N ) as a Hopf subal- 

gebra. The embedding is given by e a b | — > for any a, b = 1, . . . , N . We identify U(gl N ) 
with its image in Y(gl N ) under this embedding. It is clear from relations (4.2) that for any 
a,b=l,...,N, 

(4.6) [E ab ® l + l<g>e o6 ,T(u)] = 0. 

The evaluation homomorphism e : Y(gl N ) — > U(gl N ) is given by the rule e : T a ^ i— > e& a 

for any a, b = 1, . . . , N , and e : i— > for any s > 1 and all a, 6 . Both the automor- 
phisms p x and the homomorphism e restricted to the subalgebra U(gl N ) are the identity 
maps. 

For a gl^-module M denote by M(x) the y(gljy) -module induced from M by the 
homomorphism eo p x . The module M(x) is called an evaluation module over ^(gljy) with 
the evaluation point x . 

Denote by Y + (gl N ) the left ideal in Y(gl N ) generated by the coefficients of the series 
T ab (u) for 1 «C b < a < N. If A,B e Y(gy and A — B E Y + (gl N ) , then we will write 
A ~ S . Relations (4.2) imply the following lemma. 

Lemma 4.1. 

a) For any a, b = 1, . . . , N, the coefficients of the series [T aa {u) , T bb (y)] belong to Y + (gl N ) . 

b) For any Z G Y + (gl N ) and any a = 1, . . . , N , the coefficients of the series ZT aa (u) he- 
long to Y + (g[ N ) . 
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Therefore, for any a = 1,...,N, the coefficients of the series T aa (u) act on the space 
Y(gl N )/Y + (gl N ) by multiplication from the right, and those actions commute. 

Lemma 4.2. The coefficients of the series A(T aa (u)) — T aa (u) <g> T aa (u) , a = 1, . . . , N , 
belong to Y(gl N ) <g> Y+(gl N ) + Y + (gl N ) Y(gl N ) . 

The proof is straightforward, cf. (4.4) . 

A vector v in a Y(gl N ) -module is called singular with respect to the action of Y(gl N ) 
if Y + (gl N ) v = . A singular vector v that is an eigenvector for the action of T 11 (-u), . . . , 
Tnn(u) is called a weight singular vector. 

Example. Let M be a gl^-module and v G M a singular vector of weight (A 1 , . . . , A^) . 
Then for any x G C , the vector v is a weight singular vector with respect to the action of 
Y(gl N ) in the evaluation module M(x) . 

There is a Hopf algebra anti-involution w : Y(gl N ) — > Y(gl N ) that is defined in terms of 
the series T ab (u) by the rule 

(4.7) uj(T ab (u)) = T ba (u), a,b=l,...,N. 
In particular, for any X G Y(gl N ) we have 

(4.8) A(z*7p0) = ((w®w)(A(X))) {21) . 

Consider the series T (fc,fc+i)(V) . . . T^' k+1 \u -k+1) with coefficients in End(F® fc ) <g> 
y(fl[jy) . Relations (4.3) and Lemma 3.1 imply that 

(4.9) A^T fc) T^ k+1 \u -k + 1)... T^ k+1 \u) = 



= T<*'* +1 >(«) . . . T^ k+1 \u -k + 1) Af k} - k) . 



Therefore, the coefficients of the series T^ k ' k+1 \u) . . . T( 1 ' fc+1 ^('u — k+1) preserve V Ak . 
Define the series T Ak (u) with coefficients in End(V Ak ) ® Y(gl N ) : 



(4.10) T Ak (u) = T^ k ' k+1 \u) . . . T^' k+1 \u - k + 1) | (1 ; 



.k) 



In particular, T A1 (u) = T(u) . Formulae (4.3) and (3.3) imply that the introduced series 
satisfy the following relations: 

(4.11) (R Ak ' A \u-v)) {12 \T Ak (u)f 3 \T Al (v)f 3) = 

= (T Al (v)f 3) (T Ak (u)) {13) (R Ak ' Al (u-v)f 2) . 

Lemma 4.3. (id ® A) (T Ak (u)) = (T Ak (u)) {13) (T Ak (u)) {12) . 

Proof. The statement follows from the coproduct formula (4.5) and the definition (4.10) of 
T Ak (u). □ 
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Since the space V AN is one-dimensional, the coefficients of the series T AN (u) belong to 
the Yangian Y(gl N ) . The series T AN (u) is called the quantum determinant and is denoted 
by qdet T(u) . Explicit formulae for the quantum determinant that follow from (4.9) and 
(4.10) are 

(4.12) qdet7» = sign(a) ^ sign(r) T auTl (u) . . . T a „ tTN (u - N + 1) 

tES n 

= sign(a) si S n ( r ) T r N ,a N {u - N + 1) . . . T Tu(Tl {u) , 
beS N 

where the permutation a G Sn is arbitrary. Relations (4.11) and Lemma 3.2 imply that 
the coefficients of qdetT(u) are central in Y(gl N ) . 

Theorem 4.4. [MNO] The coefficients of qdetT(it) are free generators of the center of 
the Yangian Y(qI n ) . 

Let tr w : End(W) — > C be the trace map. For any Q e End(F) define the series 

(4.13) 7 k , Q (u) = (tr yAk ® id) (Q Ak T Ak (u)) , 

k= l,...,iV, with coefficients in Y(gl N ). We set 7 ^q(u) = 1 by definition. It is easy 
to see that 7n j q(u) = det Q qdetT(it). The series 7i t Q(u), . . . , 7n,q(u) are called the 
transfer matrices. They have been introduced in [KS]. 

Proposition 4.5. [KS] We have [7 kjQ (u) , 7 ljQ (v)] = . 

Proof. The statement follows from formulae (4.11) and (3.4): 

7 k , Q {u)7 hQ {v) = (tr yAfe ®tr yAi ®id)((g Afe T Afc («)) (13) (g A 'T A/ (,;)) (23) ) = 

= (tr VAfc ®tr yAl ®id)((Q^ = 

= (tr^^tr^^id)^^)^^^^!^'^^-^)^)" 1 X 

x (T Ak (v)f 3 \T Al (u)f 3) (R Ak > Al (u-v)) {12) ) = 

= (tr yAfc ®tr yA! <g,id)(^ x 

x (T Ak (v)f 3) (T Al (u)f 3) (R Ak ' Al (u-v)) {12) ^ = 

= (tr yAfc ®tr yAi ®id)((Q Afe )^ = 

= (tr yAfe ®tr^®id)((Q^T^(z;)) (23) (Q^T Afc ( W )) (13) ) = 

= 7i tQ (v)7 ktQ (u) . 
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Theorem 4.6. [NO] If the matrix Q £ End(V) Las a simple spectrum, then the coeffi- 
cients of the series 7i^q{u), . . . , 7n-i,q(u) and qdet T(-u) generate a maximal commutative 
subalgebra in Y(gl N ) . 

Recall that the Yangian Y(gl N ) contains U(gl N ) as a subalgebra, the embedding being 
given by e ab i-> . 

Proposition 4.7. [KS] Let Q be the identity matrix. Then for any k = 1,...,N, the 
coefficients of 7k,Q(u) commute with the subalgebra U(gl N ) C Y(gl N ) . 

Proof. Similarly to (4.6) , for any a, b = 1, . . . , AT, we have 

[eab\ v ^ k ® l + l®e ab ,T Ak (u)} = 0, 
which implies that for the identity matrix Q , 

[e ab ,7 k , Q (u)} = -(tr VAfe ®id)([e ab | i/Afe ®l,T Afc (w)]) = 

due to the cyclic property of the trace. □ 
Lemma 4.8. For any m = k, . . . , N , and any distinct i±, . . . , ik G { 1, . . . , m} , we have 

= (V ® id)(Q (<l) - • • Q ilk) T^ m+1 \u) . . . T^^\u -k+1) Af m} ^) . 

Proof. If m = k and ij = j for all j = 1, . . . , k , the claim is equivalent to formula (4.13) . In 
general, we are also using the relations 'T^ m+1 \u) = T^ m+1 \u)Q^P^ and 

P( l ^A{ m y = —A{ m y = A{ m yP( l rt , 1 ^ i < j ^ m , for the flip maps, and the formula 

(id® fc ® tr yMm _ k) ) A {m} = kl K {N _ k) \ A {k} . D 

Consider the algebra Y(gl N ) d u ]] of formal differential operators; the indeterminates 
u~ x and d u do not commute obeying the following relation d u -u~ s = u~ s d u — s-u _s_1 
instead. Following [T] , for any m = 1, . . . , N, set 

(4.14) ® m , Q (u,d u ) = 

= (tr^ m ® id) ((1 - Q^T^ m+1 \u) e~ d ") . . . (l - Q^T^ m+1 \u) e~ 9 -) 
and D ,q(w,9 u ) = 1 . 

Proposition 4.9. [T] For any m = 1, . . . , N, we have 

(4. 15 ) W „A) = g^gzgj ,•„(„) e -^. 
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Proof. Expand the product under the trace in formula (4.14), 



3>m,Q(«A) = E E 

k=0 l^ii<...<i fc ^m 



fe x 



x (tv v ® m <g> id)^* 1 ^* 1 -" 1 * 1 ^) . . . gfe)T (lfe ' m+1 )(w - jfe + 1) Af m} m) ) e" fc ^ 

Then Lemma 4.8 yields the claim. □ 

Remark. Formula (4.15) for m = N shows that the expression ® N,Q(u,d u ) serves as a 
generating function of the transfer matrices 7i,q(u), . . . , 7n,q(u) '■ 

N 

(4.16) 2W«A) = E(- 1 ) fcT ^W e " fc9u - 

fc=0 

Remark. Taking into account formula (4.9) and the equality (A^ m y) 2 = A^ m y , one can show 
that 

(4.17) (1 - Q^T^ m+1 \u) e" 9 ") . . . (l - g("») T (m,m+i) ^ e -s„j A^ m} m) = 

= Af m} m) (1 - Q( 1 'r( 1 ' m+1 '( U ) e" 9 ") ... (1 - Q( ro ) r (m, m +i) (u) e -3„) A^ m} m) . 

Since A^} is a one-dimensional projector, formula (4.17) implies the following proposition. 
Proposition 4.10. We have 

(id^SW.QK^)) = 

= (1 - Q^T^ N+1 \u) e" 9 ») ... (1 - Q^T^ N+1 \u) e~ d ") Af N} N) . 

Let be the matrix transposed to Q ; for example, E^ ah = E(, a . 
Proposition 4.11. We have w(7k,Q(u)) = 7 kt Qt(u) , where w is given by (4.7). 
Proof. The claim follows from Lemma 4.8 for 7k,Q(u) with m = k and ij = k — j + 1 , 
formula (4.9), two more observations, w(T(u)) = (T(-u))^ and (A^'j'^)^ = A^'j' fe \ and 
the standard properties of the trace: 

m(7k,Q(u)) = 

= ^ (tr^ ® u,)(Q™. . . Q^T^Xu) . . . T^\u -k + 1) A<i>*>) = 
= ^ (tr^ ® W ) (Q<*>. • • Af fc} - fc) T^ 1 ) (« - k + 1) . . . T(*'*+ 1 )(«)) = 
= ^ (tr^ ® id) (Q(*>. . . A&y *> (T( fc - fe + 1 )(«)) t . . . (T^+'Xu - k + l)) f ) = 

ll 



= I (tr^ ® id) ((Q<*>. • • A&y *> (T^^ix))* . . . (T( 1 '*+ 1 )(« - k + = 
= ^ (tr^ ® id) (T(*-*+ 1 )(«) . . . T^ 1 ' fe+1 ^('u -k+1) Afe* (Q^f . . . (<5 (1) ) f ) = 
= ^ (tr^ ® id) ( (Q«) f . . . (<5 (1) ) f T<*' fc +i) (u) . . . T^ 1,k+l \u -k + 1) A{ij*>) = 

5. The Bethe ansatz 

In this section we formulate the main technical result of this paper, Theorem 5.4. Equality 
(5.8) is, in a sense, a universal formulation of the eigenvalue problem for Bethe eigenvectors 
of the transfermatrix 7k ) Q{u) . Given a Y(gl N ) -module M and a weight singular vector 
v G M , we will apply formula (5.8) to v and obtain the Bethe eigenvector for the action 
of 7k,cj{u) in M and the formula for the corresponding eigenvalue. In Section 6 we will do 
it in detail for M being a tensor product of highest weight evaluation modules and v G M 
being the tensor product of highest weight vectors. 

Let £ = . . . , £, N ~ l ) be a collection of nonnegative integers. Set £ <a = + . . . + 
a = 1, . . . , N , and |£| = + . . . + £ Ar_1 = £ <Ar . Consider a series in |£| variables t}, ... , 

t|i , ... , if -1 , • • • , with coefficients in Y(gl N ) : 



(5.i) i^,...,^) = (tr^^, <8> id) (r^'i^i+Dcti) . . . x 

x n * (€<b+i '* <0 ^ 

(<M)<(6,j) 

where the pairs are ordered lexicographically, (a,i) < (b,j) if a < 6, or a = 6 and i < j , 
and the product is taken over all two element subsets of the set {(c,k) \ c = 1, . . . , N — 1, 
fe = 1, . . . , £ c } . Further on, we will abbreviate, t = (t\, . . . , t^v-l) . 

- N 
Lemma 5.1. Tie series B € (t) is divisible by U U (t$-t?+l) in 

Y(gl N ) [tl (t^- 1 ,)" 1 ]] • ^ 1 * < ' <€B 

Proof. Using the Yang-Baxter equation (3.2) we can rearrange the factors in the product of 
.R-matrices in formula (5.1) such that the product f| j [ R^ <a+: >'^ <a+l \tj — tf) 

appears on the right. Then formula R(u) E bc ® E bc = (u + 1) E hc ® E bc yields the claim. □ 
Let 9 : Y(gl N ) — > Y(gl N ) /Y + (gl N ) be the canonical projection. 

^ N-2 N-l C 

Lemma 5.2. The series 0(%(t)) is divisible by U U II fl (*? ~ *?) in 

a=l b=a+2 i=l j = l 

(Y(gl N )/Y + (gl N )) [t\,..., t» N z\] [[(t})" 1 , . . . , (f^)" 1 ]] . 
The proof is given in Section 11 after Proposition 11.3. 
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Set at_i r c b 1 

(5.2) B £(t ) = B 4(t ) n n i^ri n nn?^. 

a=l l^i<j^S, a 3 1 l^a<b<N i=l j = l 3 1 

cf. (5.1). To indicate the dependence on N, if necessary, we will write (£) . 

Remark. The expression B^(t) serves as a creation operator to the universal weight function 
and Bethe vectors in tensor products of evaluation modules. See more details in Section 6. 

Example. Let N = 2 and f = . Then B € (t) = Ti 2 (^) . . .Ti 2 (tji) . 
Example. Let iV = 3 and f = (1, 1) . Then 

B c (f) = T 12 (t})T 23 (t 2 ) + T 13 {t\)T 22 {t\) . 

Example. Let N = A and £ = (1, 1, 1) . Then 

B € (*) = T 12 (^)T 23 (t 2 )T 34 (t?) + 

+ ^- T T 13 (^)T 22 (t?)T 34 (t?) + 7^T l2 (t\)T 2A (tl)T 33 (tl) + 
+ (tg.^.^i) (ri 4 (tl)T 22 (^)T 33 (tf) +T 13 (t})T 24 (^)T 32 (tf)) + 



In general, B^(t) is a sum of the products 

(5.3) ^oi,6i(*i) • • • T a]6]j b ]6] (t^v~\) p(t\, . . . , t^~_\) n nfi A 

lsCa<b<7V i=l j = l 3 i 

with various a±, . . . , a^| , &i, . . . , 6|^| , and polynomials p(t\, . . . , t^lr" 1 !) . 

The direct product of the symmetric groups S^i x . . . x S^n-i acts on expressions in |£| 
variables, permuting the variables with the same superscript: 

(5.4) a l x...xa N - 1 :f(ti... i e N -_\) ^ /(£,...,& ; ... ; t"-_\, . . . , ) , 
where a a G S^a , a = 1, . . . , N — 1 . 

Lemma 5.3. [TV1, Theorem 3.3.4] The expression B^(t) is invariant under the action 

(5.4) of tie group S^i x ... x S^jv-i . 

JV 

Fix a diagonal matrix Q = ^ Q a E aa ■ Introduce the series 

a=l 

(5.5) j| : yt) = Q a T aa ^) n n^-^-i) n w -*r *) + 

i=i i=i j'=i 

+ ga+iT 0+1 , a+1 (t?) n w - o n w - ^ a + !) n ^ - ^ a+i - *) > 

j=i j =i j =i 
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a = 1, . . . , N - 1 , i = 1, . . . , £ a . Here we set £° = f N = . 

Given the following data: integers a\, . . . , a^| + / c _i, bi, . . . , b\^\ + k-i , c e {1, . . . , N} and 
i e {1, . . . , £ c } , a sequence si, . . . , + , which is a permutation of the sequence u, . . . , 
u — k + 1, t\, . . . , tf, . . . , tf^jv- 1 ! , where the variable t\ is omitted, and a polynomial t{, 
. . . , ) , consider the product 

(5.6) T ai)6l (si) •••r 0|€|+fc _ 1> 6 |4|+fc _ 1 (s| € | +fc _i) x 



- p(«;«i ^--'o n n (^f n ^) nnn^ 



a\2 



fa 11 +a_+a/llllll (+ a + i _ +a\ 
3 l^j<Kt° 3 l/ o=l j=l 1=1 ^ V 



where the factors (it — t°j) 1 in (5.6) have to be expanded as power series in u 1 . We 
denote by I^k,Q the span over C of all products (5.6) with various ai, . . . , a\^\ + k-i , bi, 
. ..,&|£|+fc-i, c,i, si, . . . ,S| € | +fc _i , and p(u; t{, . . . , t^ N ~_\) . 
Set 

£ a a— 1 £ a 

(5.7) X£ g («;f) = Q a T a » J] ~ " - \Ji ~ U , 

1=1 1 1=1 1 

a = l,...,iV. Expressions (5.7) are regarded as power series in w -1 with coefficients in 

Y( a i N )[tl...,t£ N z\]. 

Theorem 5.4. Let Q he a diagonal matrix. For any k = 1, . . . , AT, we have 

(5.8) T fc>g («)B € (t) ~B € (f) ^n X ?;Q( w - r + 1 ^) +1X C^,Q( w ^)' 

a r=l 

where the sum is taken over all k-tuples a = (ai, . . . , afc) such that 1 ^ a± < . . . < ak ^ N , 
and 11^ k q(u;t) belongs to I^,k,Q , cf. (5.6). By Lemma 4.1, the order of the product of 
elements X|g is irrelevant. 

The theorem is proved in Section 11. 

k 

Remark. Similarly to formula (4.16), the expressions Yl Yl ^£ t q( u — r + l;t) can be ob- 

a r=l 

tained as coefficients of a suitable formal differential operator: 

(5.9) Wlt, Q (u,d u ) = (l-Xl Q (u;t)e- d -)...(l-X£ Q (u;t)e- d ") = 

N k 

= £("!)* E nXfofc-r + lsOe-***. 

k=0 l^a 1 <...<a k ^N r=l 

Remark. Since the coefficients of the series qdetT(u) are central in Y(gl N ) , the first for- 
mula of (4.12) for the identity permutation a implies that 

N 

qdetT(w) B c (t) ~ B € (t) JJ T aa (u - a + 1) . 

a=l 

Together with the relation 7n,q(u) — det Q qdetT(u) , this proves Theorem 5.4 for k = N. 
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6. Tensor products of evaluation modules 

In this section we consider the action of the transfer matrices Ti 5 q( , u), . . . , Tjv,q(w) in a 
tensor product of evaluation modules over the Yangian Y(gl N ) . We will apply Theorem 5.4 
to get eigenvectors and eigenvalues of the transfer matrices. In Section 7 we will also show 
that the operators of the action are symmetric with respect to a certain symmetric bilinear 
form, see Theorem 7.6. 

Let Mi, . . . , M n be gl^-modules. Consider the tensor product M 1 (zi) ® . . . <g> M n (z n ) 
of evaluation Ffg^-modules. The operators T a b(u)\ , a,b=l,...,N, 

\M 1 (z 1 )®...®M n (z n ) 

are End (Mi ® . . . <S> M n )-valued rational functions in u, zi, . . . , z n with the denominator 

n 

II («-*), and T ab (u)\ Mi{zi) ^ MnM =S ab +0(u-') as u ^ oo . 

Further on, we will abbreviate, z = (z±, ... , z n ) . For any k = 1, . . . , n , the operator 

(6.1) 3JS*-*^(.;,) = W«)| M , ( „ )8 ... SM „ (2n) 

n fe— 1 

is a rational function in ix, zi, . . . , z n with the denominator ] j \\ (u — j — Zi) , and 

i=l j=0 

t m^...®m„ (w; ^ = tr ^ fc QAk + ( u -i) as u _ oo . Th e operators Tjf£®"® M »(u; z), . . . , 
c yM 1 ®...®M n z ^ are ca jj e( j ^he transfer matrices of the XXX -type model associated with 
the Lie algebra gl N . 

Let Mi, . . . , M n be highest weight g [^-modules with highest weights Ai, . . . , A n , Aj = 
(Aj, . . . , Af) , and highest weight vectors v±, . . . ,v n . Then for any z\, . . . , z n , the vector 
v\ ® . . . <g) v n is a weight singular vector with respect to the action of Y(gl N ) in the module 
M x {zi) ® . . . ® M n (z n ) and 



n 



I r 1i — Zj ~\~ A.^ 

T^M^®...®^ = ui (8) . . . (8) u n , a = 1,...,AT. 

I I w - ^ 

Let £ = . . . , ^ Ar_1 ) be a collection of nonnegative integers. The function 

jv-i C n n-2 r C +1 

(6.2) z) = B € (f) V1 ®...®i; n n nn^-^ nnn (*r 

a=l i=l j = l a=l i=l j = l 

is a polynomial in t\, . . . , , zi, . . . , z n , cf. (5.3) and Lemma 5.2. We call it the universal 
weight function for the XXX- type model. 

N 

Let Q = Yl QaE aa be a diagonal matrix. The system of algebraic equations 

a=l 

(6.3) Q a n w - zj + a?) n (*? - r 1 + !) n (*? - ^ a - *) n (** - ^ a+i ) = 

J=l 3=1 3=1 3=1 

3=£i 

t-a. — 1 /-a /-o+l 

n C ? ? 

= Qa+i n (*? - + A f ') n ^ - ^ a_i ) n (** - ^ + v n (** - - *) > 

J=l J=l 3=1 3=1 
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a = 1, . . . , N — 1 , % = 1, . . . , £ a , £° = £ N = , is called the Bethe ansatz equations. Say 
that a solution t = . . . ,t^~_\) of system (6.3) is off-diagonal if 7^ for any a = 1, 
...,iV-l, l^z<j'^e,and if for any a = 1, . . . , AT - 2 , i = l,...,£ a , j = 1, 

. . . , £ a+1 . For any a = 1, . . . , N, set 

<«> x, Q(u;t; , ; a, . n ^ ff ^ n ^ ■ 

i=l i=l 1 i=l 1 

Theorem 6.1. Let Mi, . . . , M n be highest weight &l N -modules with highest weights Ai, 
. . . , A n and highest weight vectors vi, . . . ,v n . Let Q be a diagonal matrix and t = (t\, 
. . . , t^r" 1 !) an off-diagonal solution of system (6.3) . Then for any k = 1, . . . , N , we have 

k 

(6.5) 7^®-® M -{u;z)W } ^-® v -{t-z) = B^ 1 ®-®*" (t; z) H %^ Q (u - r + l;t; z; A) , 

a r=l 

where the sum is taken over all k-tuples a = (ai, . . . , ajt) such that 1 ^ ai < . . . < at ^ N . 
Proof. Theorem 5.4 and formula (5.6) imply that 

^®-® M "(u;z)Bj 1 ®-® t, "(t;z) = 

k 

= M v 1 ®...®v n{t . z) J- 11 X?; g («-r+l;t;z;A) + U£®g®^ («; f; z; A) , 

a r=l 

where K^ 1 ® g ® w " (it; t ; z; A) is a ratio of a polynomial with coefficients in Mi ® . . . <g) M n 
vanishing at solutions of system (6.3) and the product 

N-l N-2 C C +1 

n n n nnr-^) 

a=l l^i<j^£ a a=l i=l J = l 

that does not vanish at off-diagonal points (t\, . . . ,t^ N ~_\) . This yields formula (6.5). □ 

For an off- diagonal solution t = (t{, . . . , t* N ~_\) of system (6.3) , the vector B^ 1 ®-®^™ (t; z) 
is called the Bethe vector. 

k 

Remark. Similarly to formula (5.9), the eigenvalues Yl II q(w — r + z; A) of the 

transfer matrices 7 k g" n (u;z) can be obtained as coefficients of a suitable difference 
operator: 

(6.6) ^, Q (u,d u ;t;z;A) = (l - %l Q (u; t; z; A) e~ 9 -) . . . (l - X% Q (u; t; z; A) e~ 9 -) = 

N k 

= E(- 1 )" E Xm; Q {u-r + l;t;z;k)e- kd ". 



k=Q l^a 1 <...<a k ^N r=l 
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Remark. Since the coefficients of the series qdetT(it) are central in Y(qI n ) , the first for- 
mula of (4.12) for the identity permutation o implies that 



n N ~ . „ 

^etT(u)\ , ^ , = TT TT u - a+1 -« + A * 

M v ' \M-, z-,)(^...^M„ z„) 11 11 



lMi(zi)®...®M„(«„) 1111 u-a+l-z 
i—1 a=l 

n N 

Together with relation 7n : q(u) = det Q qdetT(u) \ j [\ (u — a + 1 — z^) , this proves 
Theorem 6.1 for k = N. 4=1 a=1 

Proposition 6.2. [KR1] Let Q be the identity matrix, that is, Q\ = . . . = Qn = 1 , and 

t = it\, . . . ,1^^) an off-diagonal solution of the Bethe ansatz equations (6.3). Then the 
Bethe vector M V £®'"® Vn (t; z) is a singular vector with respect to the gl N -action in Mi® 



n n 



...®m„ of weight (eaj-c 1 , EA?+e i -e 2 ,...,EAf r +e jv - 1 ). 

i=l i=l i=l 

Proof. Clearly, for any t the vector B^ 1 ®'"®^™ (t; z) has the indicated weight. The fact that 
the Bethe vector W^-®"'® Vri {t; z) is singular is proved in Section 11. □ 

Remark. Let Mi, . . . , M n be any g [^-modules, not necessarily highest weight ones, and let 
v G Mi® ...(g) M n be any weight singular vector with respect to the action of ^(gtjy) in 
the module Mi(zi) <S> ■ ■ ■ <8> M n (z n ) . Define polynomials <?i (it), . . . , q v N {u) by the rule 

T aa (u)v Y[(u-Zi) = q v a {u)v , a = 1, . . . , N. 

i=l 

The function 

AT-l C a n iV-2 C C +1 

b%+,z) = M^t; Z )v n n Uitt-zj) nnn cr 1 -*?) 

a=l j=l i=l a=l i=l j = l 

lJV-1 



is a polynomial in t\, . . . , t ?JV _i , cf. (5.3) and Lemma 5.2. For a = 1, . . . , N, set 

xa(«;*> - «.«:(-) n ^ n" fi ^ ■ 

1=1 1=1 * 1 = 1 ' 

Similarly to Theorem 6.1, we can show that for any off- diagonal solution t of the system 

^a— 1 £a+l 

Qa «) n w - t 1 + !) n w - ^ a - !) n ^ a - o = 

j=i j =i o= i 

^a—l £a ^a+1 

= Q a+ i q v a+1 (t?) n (*? - t 1 ) n (** - ^ + !) n w - *r 1 - *) > 
j=i j=i j=i 
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a = 1, . . . , N - 1 , i = 1, . . . , £ a , where £° = £ N = , we have 

T Jg8-®Wn (ti;z)B?(f;z) = B » (f;z) £ X &>-r + l;f). 

a r=l 

Moreover, for any t , 

N-i , n i \ 

qdetT(n) B\{t;z) = B\(f,z) J [ U a {u-j) \ [ _._ . 

j=0 v i=i w J 

7. Shapovalov form 

Define an anti-involution r : C/ (fltjv) — *■ ^ (fltjv) by the rme r ( e afe) = e ba , a, 6 = 1, ... , 
iV . Let M be a highest weight g [^-module with highest weight vector v . The Shapovalov 
form Sm '■ M <g> M — > C is the unique symmetric bilinear form such that Sm(v,v) = 1 , 
and Sm(Xwi,W2) = Sm(wi,t(X)w2) for any X G U(gl N ) and any wi,u>2 G M. For an 
irreducible module M , the form Sm is nondegenerate. 

Lemma 7.1. For any X G , it?i, tU2 G M , and x G C , we nave 

The proof is straightforward. 

Let M±,...,M n be highest weight gl^r-modules. Then formula (4.8) and Lemma 7.1 
imply that for any X G Y(gl N ) and any w±, w 2 G M\ ® . . . <g> M n , we have 

(7.1) ^^..^^(Xj^^^^^^,^) = 

= s Ml ® ...®Sm„ K (-(^)| M „ (z „ ) ^...^ Ml(zi) ) (n - 1) ^) , 

where .Sm, is the Shapovalov form for the module Mj , z = 1, . . . , n . 

Proposition 7.2. For any k = 1, . . . , N, and any wi,w 2 G Mi (g) . . . <g) M n , we have 

S Ml ®...® S Mn (7^®-® M " { U ;z 1 ,...,z n )w u w 2 ) = 
= S Ml ®...®S Mn ( Wl , (T fc M Q?- 8Ml («; *„,..., zi)) ( "- 1) «;2) , 
where tie superscript f stands for the transposition of matrices. 

Proof. The claim follows from formulae (6.1), (7.1) and Proposition 4.11. □ 

We remind below some facts on /^-matrices for tensor products of evaluation modules over 
the Yangian Y(gl N ) . Propositions 7.3-7.5 follow from Drinfeld's results on the universal 
/^-matrices for Yangians [D] . 
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Proposition 7.3. Let L and M be highest weight gi N -modules with the respective 
highest weight vectors v and w. For generic x,y G C there exists a unique operator 
R LM (x -y) e End (L<g>M) such that 

(7.2) R TM (x-y)X\ = (X\ ) (21) R T M (x - y) 

V ' LMK yj \L{x)®M(y) v Im( 2 /)®L(x)'' LMV y; 

for any I £ 7(5^), and 

(7.3) R LM {u) v <S> w = v <S> w . 

The operator R LM (u) is called the rational R-matrix for the tensor product L®M . 
Example. Let L = V M and M = V Am . Then 

r lm {u) = R Ai > Am (u) ^ +max(m '^ 0) n Tr 1 • 

i=0 j=0 J 



Proposition 7.4. R LM (u) = [{R ML (-u)) {2l) ) . 



Proposition 7.5. Let Mi,M 2 ,M 3 be highest weight gl N -modules. Then the operators 
R M , M ,{u) , 1 ^ % < j ^ 3 , satisfy the Yang-Baxter equation 

d( 12 ) i \ d( 13 ) ( \ d( 23 ) i \ d( 23 ) i \ d( 13 ) / ^ d( 12 ) i \ 

K M 1 M 2 y u i - u v k m 1 m 3 v u v k m 2 m 3 k u v - K M 2 M 3 y u 2) K M 1 M 3 y u i) K M 1 M 2 y u i - u v ■ 

By the definition of evaluation modules over Y(gl N ) , formula (7.2) means that the op- 
erator R LM (u) commutes with the {([^-action in L ®M and for any a, b = 1, . . . , N, 

( N \ ( N \ 

r lm( u ) ue ab <S> 1 + E e ac <S> e cb ) = ue ab ® 1 + J2 ) r lm( u )- 

V c=l 7 V c=l 7 

The operator R LM (u) is uniquely determined by these properties together with the nor- 
malization condition (7.3), which implies that the restriction of R LM (u) to any weight 
subspace of the gl^-module L (g>M is a rational function of u, and 

(7.4) R L m(u) = 1 + Oiu- 1 ) , u -> oo . 
In addition, for any w\ , wi e L <g) M , we have 

(7.5) 5x<g> Sm(R lm ( u ) w 1i w 2) = S L <g) S M (w 1 ,R LM (u)w 2 ) . 
Let Mi, . . . , M n be highest weight gl^-modules. Set 

(7-6) R Mi ...mS z )= II ^m?m,(^-%) eEnd(Mi®...®M n ). 
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By Proposition 7.4 we have 



(-^Mi...M„( 2 )) — ]^[ ^M j M l ( z j Zi ) ■ 

For any wi,w 2 G Mi ® . . . <g) M n , formula (7.5) and Proposition 7.5 yield 

(7.7) S Ml ® .. S Mn {R Ml _„ Mri (z)w 1 ,w 2 ) = S M 1 ®---®S Mn {w ll R Mi _ Mn (z)w 2 ), 
and for any X G Y(gl N ) , formula (7.2) implies 

(7.8) i2 M M (z) XI = (X\ ) {n " A) R M M (z). 

V ) M 1 ...M n \ ) \M 1 (z 1 )®...®M n (z n ) V \M n (z n )®...®M 1 (z 1 )' Mi...M n V 7 

Define a bilinear form Sj^ <g, (g> M on Mi ® . . . <g) M n by the rule 

( 7 - 9 ) Sm!®...^^!)^) = (S Ml ® • • • ® 5jifJ(wi, %...jif„ • 

For irreducible modules Mi, . . . , M n and generic z , the form ®M n ls nondegenerate. 

By equality (7.7), the form ^ i8 ^ Mn is symmetric. 

Theorem 7.6. For a symmetric matrix Q, Q = , the operators 7j^^"'^ Mn (u; z) , 
are symmetric with respect to the symmetric bilinear form ®M n ■ 

Proof. Proposition 7.2 and formula (7.8) yield 

OZ (<jM 1 ®...®M n( \ \ _ qz ( <T M 1 ®...®M n( S \ 

for any k = 1, . . . , N . This proves the statement. □ 

8. Current algebra fltjv[x] 

Let £jl/v[#] be the Lie algebra of polynomials with coefficients in jjtjy with the pointwise 
commutator. We call it the current algebra. In this section we construct, following [T], 
analogues of the transfer matrices 7i : q(u), . . . : 7n,q(u) for the algebra U(qI n [x]) , and 
formulate the counterparts of Theorems 5.4 and 6.1. 

We will be using the standard generators of fl(jv[x], e)2 = e a &x s , a, b = 1,...,N, 
and s = 0,1, . . . , which commute as follows: [e^ h , ] = ^&c e ad ~~ ^ade^ +s ^. We will 
identify the Lie algebra gl N with the subalgebra in gljvf^] of constant polynomials, thus 
making no difference between the elements and the generators e a & of q\ n . 

Let n + = Ce a b and n_ = Ce^a be the standard nilpotent subalgebras in 

l^a<b^N l^a<b^N 

qI n , and let n±[x] be the corresponding subalgebras in gljvf^] • If A,B G C^Cflljvt^]) an d 
A — Be U(gl N [x])n + [x] , then we will write A ~ £? . 

A vector w in a ^[^[xj-module is called singular with respect to the action of if 
n + [x]u = . A singular vector i> that is an eigenvector for the action of e\?} , a = 1, . . . , 
N , s = 0, 1, . . . , is called a weight singular vector. 
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Let M be a g [^-module. For any y G C denote by M(y) the corresponding evaluation 
gl^ [a] -module, that is, 4& | M(v ) = 2/ s e a 6| M , a, b = 1, . . . , N, s = 0, 1, . . . . 

We organize elements e)2 , a,b = 1, . . . , N , s = 0, 1, . . . , into generating series 



oo 

(8.1) L ab (w) = 5>£ }u " a-1 > a,b=l,...,N. 

s=0 

N 

and combine those series together into a series L(u) = E ab <E> L ab (u) with coefficients 

a, 6=1 

in End(C Ar ) <8> {/(fl^x]) • Notice the flip of subscripts of the generators in the sum (8.1); 
this is done to have better notational correspondence further with the Yangian case. The 
commutation relations in U(gl N [x]) can be written in the following form: 

(it - v) [L ab (u) , L cd (v)] = 5 bc (L ad (u) - L ad (v)) - 5 ad (L cb (u) - L cb (v)) , 

that amounts to the equality for series with coefficients in End (C^) <g> End (C^) <g> U(gl N [x]) 

(8.2) (u - v) [L^(u) , L< 23 >(t/)] = [L^(u) + L«*\v) , P^ 12 )] . 

Here P G End(C JV ) <g> End(C JV ) is the flip map. 

Fix a matrix K G End(y) and consider the formal differential operator 

(8.3) V K (u,d u ) = 

= {tr v9K ®M)((d u -KW-LV> N + 1 ^ A^f . 

of order N. Let Go,k(u) > • • • > Gn,k(u) G ^(flW^DH 1 ' -1 ]] be its coefficients: 

N 

(8.4) V K (u,d u ) = £(-1)* 0""*- 

k=0 

The series Gi,k(u), ■ ■ ■ , Gn,k(u) are called the Gaudin transfer matrices. 

N 

Example. Clearly, Go,k(u) = 1 because tr v<S)N A{ N y = 1. Let K = K ab E ab . 
We also have that a ' 6=1 

N 



(8.5) Q ljK (u) = (tr v ® id) (KW + L( 12 )(«)) = ]T(K aa + L a »), 

o=l 

= (tr vm ® id) (((KWu^Mj^ + L^^-L^MjAS^) = 

1 ^ 
= g ((^i,«r(«)) 2 -(^-l)^i,if(«) " (K ab + L ab (u))(K ba + L ba (u))^) 

a, 6=1 

where G = dG/du . 
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The next two propositions were obtained in [T] from the Yangian case. We will give their 
proofs at the end of Section 10. 

Proposition 8.1. We have 

A? N} N) V<£ +1 \u,d u ) = (d u -KW-L^\u))...(d u -K^-L^ N ^(u)) Af N} N) . 
Proposition 8.2. We have \Gk,K{u) , Gi,k(v)~\ =0. 

Proposition 8.3. Let K be the zero matrix. Then for any k = 1, . . . , N , the coefficients 
of Qk,K{u) commute with the subalgebra U(gl N ) C U(gl N [x]) . 

Proof. Similarly to (4.6) , for any a, b = 1, . . . , N, we have [E a b ® 1 + 1 <S) e a b , L(u)~\ = . 
Hence, the claim follows from formulae (8.3) , (8.4) , and the cyclic property of the trace, 
because [E^ + . . . + E% >,A frf>] = . □ 

We extend the anti-involution r : U(gl N ) — > U(gl N ) to the algebra U(gl N [x]) by the rule 

(8.6) T(e { a s 6 } ) = e&> , a, b = 1, . . . , N, . = 0,1,.... 

The extension is consistent with the standard embedding U(gl N ) U(gl N [x]) . 

Proposition 8.4. We have r[Qk,K(u)) = Qk,Kt(u), where the superscript f stands for 
the transposition of matrices. 

The statement is proved in Section 10. 

Let £= . . . , £, N ~ 1 ) be a collection of nonnegative integers. Set |£| = £ x + . . . + ^ Ar_1 . 
Consider an expression in |£| variables t\,. . .,t^ , ... , t^ -1 , . . . , t^Z^ with coefficients 

in U(n_[x]) : 

(8.7) F^},...,^" 1 ,) = 

f > / -j^ rn ab 

= Sym tl _ tl . . . Sym^-i^. ^ ]T II ( ^TT II L ^(%+*) X 

b-2 

n 



b-2 

X 



/ c+1 - / c 



where the sum is taken over all arrays of nonnegative integers m = (m a b) i^a<b^N such 

a N 

that ^ m cb = £ a , and = ^ m rs . 

c=l b=a+l (r,s)<(a,b) 

S>C N-2 C C +1 

Lemma 8.5. The product F^t] /,\ 1 ) jj Jj ! [ i/" ' ; tf) is a power series in 

a=l i=l j=l 

(t{) , . . . , (t^vl 1 ! ) with coefficients in U(n_ [x]) . 
The proof is straightforward. 
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N 

Let K be a diagonal matrix, K = Y K a E aa . Introduce the expressions 

a=l 

(8.8) r^ K (t) = K a - K a+1 + L aa (t<}) - L a+ha+1 (t?) + 

~^ ^ ; j-a j.a— 1 ^ ^ ; +a _ -i-a ^ ] , a j-a+1 ' 

i=l ' ~ S j = l ' J j 1 '< ~ S' 

a = l,...,iV-l, i = l,...,f a . Here we set £° = £ N = 0. 
Denote L^(it) = -—L ab (u) . Given the following data: 

Uj Hi 

i) an integer m not exceeding |£| + k — 1 , integers ai, . . . , a m , b±, . . . , b m , c G { 1, . . . , iV} 
and ie{l,...,f}, 

ii) a sequence si, . . . , s m , where each sj is either u or one of the variables t\, . . . , tf, . . . , 
t^N~-i , where the variable t\ is excluded, 

iii) nonnegative integers r±, . . . , r m such that rj = if Sj ^ u , and Yl r j^k, 

iv) a polynomial p(u; t{, . . . , t^ N _ 1 ) , 
consider the product 

(8-9) L^l( Sl ) . . . L[ r ±(s m ) dlfAt) x 

x P («i *f»- 1 ,) n n n n n 

o=l j = l 1 J J o=l j=l 1=1 *i _t J 

where the factors (u — t'j)~ k in (8.9) have to be expanded as power series in u~ x . We denote 
by J^,k,K the span over C of all products (8.9) with various m , ai, . . . , a m , bi, . . . , b m , 
fi-, • • • ■> fm 1 c,i , si, . . . , s m , and p(u; t\, . . . , t^ N _ 1 ) . 
Consider the formal differential operator 

Mt, K {u,d u ;t) = fl (d u -K a -L aa (u)-J2^^l + £ rz^r) • 

Ka^7V V i=l i=l U ^ 

Let Z^i^xiu; t), . . . , Z^n,k(u; t) be its coefficients: 

TV 

M € ,x(u,0 u ;f) = ^(-l^^MR^. 



fc=0 

Expressions ^^^(m; £), . . . , Z^n,k{u; t) are regarded as power series in u~ l with coeffi- 
cients in U(gl N [x])[t\, . . . ,t^ N ~_\] . 

Theorem 8.6. Let K be a diagonal matrix. For any k = 1, . . . , N , we have 

(8.10) Pfc,jr(«)F { (t) ~ F ? (t)Z ? , fcjK ( W ;t) + W C fc)Ar (u;t), 

where "W^ fe ^(it; t) belongs to J^k,K ■ 
The theorem is proved in Section 12. 
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Remark. It is easy to see that the coefficients of the series Gi,k(u>) are central in U (qIn[x]) 
and Z^i iK (u;t) = Gi,k{u) for any £. This implies Theorem 8.6 for k = 1 . 

9. Tensor products of evaluation modules of the current algebra 

Let Mi, . . . , M n be gl^-modules. Consider the tensor product Mi(zi) <S> ■ ■ ■ <S> M n (z n ) 
of evaluation gl^y [a;] -modules. Then 

n (i) 

(9.1) L ab (u)\ , , =y^3«_ a = l,...,N, 

i=l 

where we consider = 1®(* _1 ) <g> e a & ® as an operator acting in Mi <g> . . . <g> M n , 

cf. (8.1) . By formulae (8.3) , (8.4) , for any k = 1, . . . , n , the operator 



(9-2) G™f-® M "{u;z) = Q k ,K(u)\ Ml 



(zi)®...®M„(z„) 



n 

is a rational function in it, zi, . . . , z n with the denominator j ] (u — Zi) h and 



i=l 



G™f-® M -(u;z) = tvK Ak + 0(u~ l ), 



u — ■> oo . 



The operators ^®-® M "(u; z), . . . -® M "(u; z) are called the trans/er matrices of 

the Gaudin model associated with the Lie algebra gljy . 

Let r be the anti-involution (8.6) on U(gl N [x}) . Recall that for a highest weight gl N 
module M we denote by Sm the Shapovalov form on M . 

Let Mi,...,M n be highest weight gl^-modules. For any zi,...,z n , X e U(gl N [x]) , 
and wi,W2 G Mi ® . . . <g) M n , we have 

S Ml ® ■ ■ ■ ® S M „ (^| Ml(2l)8 ... 8M „ ( ^^1^2) = 
= S> Ml ® • • • ® S M „ K^W| MiW ,., MB(2n) ^) • 

Then Proposition 8.4 and formula (9.2) yield that for any k = 1, . . . , N, 

S Ml ®...®S Mn ( gffi-Wn {u . z) Wl , W2 ) = s Ml ® . . . ® S< M „ ( Wl , gffi- 9 "* ( u - Z ) W2 ) . 

In particular, we have the following result. 

Theorem 9.1. For a symmetric matrix K, K = , the operators Gff^'"® Mn (u; z) , 
k = 1, . . . , N, are symmetric with respect to the tensor product Smi® ■ ■ ■ <S> Sm„ of the 
Shapovalov forms. 

Let Mi, . . . , M n be highest weight g [^-modules with highest weights Ai, . . . , A n , Aj = 
(Aj, . . . , Af ) , and highest weight vectors . . . , i; n . Then for any zi, . . . , z n , the vector 
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v\ <S> ■ ■ ■ <S> v n is a weight singular vector with respect to the action of fltjvl^] in the module 
Mi(zi) <g> ...®M n {z n ) and 

A<? 

Laa{u) v\ <S> ■ ■ ■ <8> v n = vi <g> ■ ■ ■ <E> v n — , a = 1, . . . , N. 

U — Z; 
i=l 1 

N 

Fix a diagonal matrix K = K a E aa . Let £ = . . . ,^ Ar_1 ) be a collection of non- 

a=l 

negative integers. Set £° = £ N = . The system of equations 

(9 .3) y ! h ! h + y — 1 _ 2 y -J— + y — 3_ = K a+1 -K a , 

j 1 ' 'I j=l ~ '.i j=l / j=l '> ~ '.i 

3=£i 



a = 1, . . . , N — 1 , i = 1, . . . , £ a , is called the Bethe ansatz equations. We always assume 
that for a solution t = (t\, . . . , t^X) of system (9.3) any denominator in equations does 
not equal zero unless the corresponding numerator equals zero. 
Consider the differential operator 



(-a — 1 /-a 

— > / " »« « i « i \ 

(9.4) M„(M.^;A) = J] (^-^-E^TF "Et-^i +E^j 

I^aA i=x U Z * i =i U ~ t i i=i 

Let Z^i^k{u] t; z; A), . . . , Z^n,k(w, t; z; A) be its coefficients: 

N 

(9.5) M^(«,9 u ;^;z;A) = ^ (-l) fc Z Ck , K (u; t; z; A) a 



AT 

fc=0 



The expression (8.7) defines a vector-valued rational function ]pj 1 ®---® v «(^ z } Q f the vari- 
ables t]_, . . . , tfjv-i 5 ^1) • • • j ? 



(9.6) F^-^"(t;^) = F € (t) Vl (8)...(8)i; n , 

with the denominator J] f] f] (£? - ^) f] f] J] (t" +1 - ^ a ) , cf. Lemma 8.5. The 

i=l a=l j = l a=l i=l j=l 

function F^ 1 ®'"®'"™ (t; 2) coincides with the universal weight function used in integral formu- 
lae for hypergeometric solutions of the Knizhnik-Zamolodchikov equations, see [Ma] , [SV] , 
[RSV]. 

Theorem 9.2. Let M±, . . . , M n be highest weight gi N -modules with highest weights Ai, 
. . . , A n and highest weight vectors v±, . . . , v n . Let K be a diagonal matrix and t = (t\, 
. . . , t^l\) a solution of system (9.3). Then for any k = 1, . . . , N, we have 

(9.7) gMf---® M -(u;z)F v ^---® v "(t;z) = Z^ k , K (u; t; z; A) (t; z) . 
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Proof. The statement follows from Theorem 8.6 and formulae (9.2), (9.6). □ 

Remark. Since the coefficients of the series Gi,k(u) are central in U(gl N [x]) , formulae (8.5) 
and (9.4) yield 

a=l ^ i= i U Zi J 

which implies Theorem 8.6 for k = 1 . 

Proposition 9.3. [RV] Let K be the zero matrix, that is, K\ = . . . = Kn = , and t = 
(£[,..., t^X) an off-diagonal solution of the Bethe ansatz equations (9.3) . Then the Bethe 
vector Y V £®'"® Vn (t; z) is a singular vector with respect to the gl N -action in Mi <g> . . . ® M n 



n 



of weight (EAJ-C 1 , EAf+^-^-.^EAf+r- 1 )- 

i=l i=l i=l 

Proof. The claim follows from identites for the universal weight function [RSV] and quasi- 
classical asymptotics of hypergeometric solutions of the Knizhnik-Zamolodchikov equations 
[RV]. " □ 

Remark. Let Mi, . . . , M n be any g [^-modules, not necessarily highest weight ones, and let 
v G Mi® . . . <g) M n be any weight singular vector with respect to the action of fll/vM i n 
the module Mi(zi) ® . . . <g> M n (z n ) . Define rational functions r\(u), . . . , by the rule 

Laa{u) v = r%(u) v , a = 1, . . . , N , and consider the differential operator 

Ml K (u,d u ;t) = J! (du-K a -rl(u)-J2-^l + £ ^) • 

i^scaA i=i i=i w r ^ 

Let Z| 1 K (u;t; z), . . . , Z£ N K {u;t; z) be its coefficients: 

N 

Ml K (u,d u ;t;z;A) = £ (-l) fc Z^u; t; z) d?~ k . 

k=0 

Similarly to Theorem 9.2, for any solution t of the system 

i ^ - i ^ ^ 

r a(^) + ^2 7^ I^T ~~ 2 ^2 +a _ +a + ^2 Ta I^+T = Ka+1 ~ Ka ' 
3 = 1 ~ ' l j=l ' ' / j=l 'i ~ '.i 

a = 1, . . . , N — 1 , i = 1, . . . , £ a , we have 

QMf...®M n{u . z)¥i{i)v = ZlKK{u .~ t . z ) ¥i (t)v. 

Moreover, for any £ and t , 

N 



a =i 
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10. Filtration on Y(gi N ) 

In this section we are going to relate the series Ti^u), . . . , 7n,q(u) with coefficients in 
the Yangian Y(gl N ) and the series Gi,k(u), . . . ,Gn,k(u) with coefficients in the algebra 
U(gl N [x]) , and to prove Propositions 8.1-8.4. The results of this section are essentially 
borrowed from [T] . 

Given a filtered algebra A with an ascending filtration . . . C A s -i C A s C . . . C A, 
possibly double infinite, we denote by grad^ : A s — > A s /A s -i the canonical projection. We 
will identify the quotient spaces A s /A s -\ with the corresponding homogeneous subspaces 
in the algebra 

grad.4 = (J) Ar/A r -i 

and will regard grad^ as a map from A s to grad A . Abusing notation, we will write just 
grad s , dropping the superscript, if it can cause no confusion. 

For a filtered algebra A , we consider the algebra End (V) <8> A to be filtered by the 
subspaces (End (V) <g> A) s = End (V) <g> A s . 

The Yangian Y(gl N ) admits a degree function such that deg = s — 1 for any a, b = 
l,...,iV, s = 1,2, see Section 1.20 in [MNO]; the function deg here coincides with 
the function deg 2 in [MNO]. Then Y(gl N ) is a filtered algebra with an ascending filtration 
Y(gl N ) C Y(gl N ) 1 C . . . C Y(gl N ) , where the subspace Y(gl N ) s contains elements of degree 
not greater than s . 

Theorem 10.1. [MNO, Theorem 1.26] We have grad(Y(gl N )) = U(gl N [x]) . 
Moreover, grad^^T^) = e£ } , a, b = 1, . . . , N, s = l,2,.... 

The anti-involution w is compatible with the filtration on Y(gl N ) . For any X e Y(gl N ) 
we have deg zu(X) = deg X , and if deg X = r , then 

(10.1) grad r (ro(X)) = r(grad r X) . 

We extend the filtration from the Yangian Y(gl N ) to the algebra Y(gl N ) [[w _1 , d u }} by the 
rule deg u~ x = deg d u = — 1 . Notice that the indeterminates u~ x and <9 U do not commute 
obeying the following relation d u -u~ s = u~ s d u — su~ s ~ l instead. Clearly, 

grad(Y( y[[ W - 1 ,^]]) = U{g[ N [x])[[u-\d u ]\. 

It is easy to see that the series T a b(u) — 5 a b £ Y(gl N )[[u~ 1 , d u ]] has degree —1 and 

(10.2) grad.! (T ab (u) - 5 ab ) = L ab (u) , a,b=l,...,N. 

Further on we assume that in the definition (4.13) of the transfer matrices Q is a series 
from End(V)[[£]] instead of being an element of End(y). Then the obtained transfer 
matrices 7^q{u), . . . , 7n,q(u) are power series in u~ l and ( with coefficients in in Y(gl N ) . 
We also regard them as elements of Y(gl N )[[u~ 1 , d u , (]] . The results of Sections 4 and 5 
naturally generalize to the described setting. We extend the filtration from Y(gl N )[[u~ 1 , d u ]] 
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to y(flljv)[[w \<9 U ,C]] by setting deg ( = — 1. Similarly, in formulae (8.3), (8.4) we will 
assume K being a series from End(V)[[£]] . 

Recall that 7 0i q(u) = 1 by convention. For any k = 0, . . . , N, set 
(10.3) = ^ ± (-!)« ^ %, Q (n) . 

For example, $o,q(u) = 1 and Si,q(w) = Ti ; q(w) — AT. Formulae (10.3) for all k = 0, . . . , 
iV taken together are equivalent to the identity 

N N 

(10-4) £("!)* = £(-1)^(^+1)^. 

fc=0 z=o 

The series §i 5 q(m), . . . , Sjv,q(w) are called the modified transfer matrices. 

Theorem 10.2. [T] Let deg(Q - 1) < -1 . Then deg (B kjQ (u)) = -k for any k = 0, . . . , 

AT. Moreover, 

(10-5) grad_ fc (S fc , g (w)) = Gk,K{u) 

where K = grad_ x (<5 — 1) ■ 

Proof. Taking the (N — m)-th derivative of formula (10.4) with respect to y and making 
there the substitution y = e du — 1 , we get 

( 10 - 6 ) E^T^fr^W^- 1 ) 



cf. (4.15). Exploiting the standard identity 

yj-iy 

^-^ rl (s — r 

r=0 V 

we obtain that 

(!0-7) B k , Q (u) = EC" 1 )™ ff^T V m , Q (u,d u )e^(e^-l) k -™. 

Since deg(T(w) - l) = -1 , deg(e a " - 1) = -1 , and deg(Q - 1) < -1 , formula (4.14) 
implies that deg(%) m ,Q(u,d u )} = —m. Hence, deg (§A;,q(w)) = —k by formula (10.7). 

Observe that grad_ 1 (T(w) - l) = L{u) , cf. (10.2), and grad.^e 9 " - 1) = d u . Then 
computing grad_ m D mj Q(w, d u ) in two ways using either formula (4.14) or formula (10.6), 



Em) 



1=0 



i (N-iy. 

(m-l)l 



7i,Q(u)e 



(m—l)d u __ 



(N-m)\® m , Q (u,d u )e md -, 



0, 
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we get 

(tr^ m ® id) ((d u - KM - L^ m+1 \u)) ... (d u - - L< m > m +V(u)) A^) = 



(iV-m)! ^ ' (m-k)\ 

The last step is to compare the obtained formula for m = N with formulae (8.3) and (8.4) , 
which completes the proof. □ 

Proof of Propositions 8.1-8.4- Let Q = 1 + (K . Propositions 4.5, 4.11 and 4.10 together 
with formulae (10.3), (10.5), (10.1) and the equality gmd _ n (®n,q('u>, c? u )) = T>£K(u,d u ) 
imply respectively Propositions 8.2, 8.4 and 8.1 with the matrix K replaced by . Since 
all the expressions Gi,c,k{u), . . . , Gn,ck(u) and ©^(u, <9 U ) are polynomials in C, we can 
evaluate them at £ = 1 , which completes the proof. □ 

11. Proof of Theorem 5.4 

We will prove Theorem 5.4 by induction with respect to N . The key points are Propo- 
sition 11.3 and 11.5, and formula (11.25). 

Set W = C^ -1 . Let wi, . . . , Wjv-i be the standard basis of W and Vi,...,Vjv the 
standard basis of V = C N . We identify W with its image in V under the embedding 
w a i-> v a+ i , a = 1, . . . , N - 1 . 

Let P^" 1 ) e End(17® 2 ) be the flip map and R^'^iu) = u + P^" 1 ) the rational 
.R-matrix used in the definition of the Yangian Y(gl N _ 1 ) . The .R-matrix R(u) preserves 
the subspace W® 2 C V® 2 and the restriction of R(u) on W® 2 coincides with . 

We will use the embedding V : ^(flU-i) ^ ^(flU) : V>( T if _1> («)) = T a+hb+1 {u) , 
a, 6 = 1, . . . ,N — 1 , where T a ^ ^(it) are series (4.1) for the algebra Y(qI n _ 1 ) . It is clear 

that vft(flU-i)) c . 

Let W(x) be the evaluation y(g[ Ar _ 1 )-module and 7r(x) : ^(fltjv-i) — End(VF) be the 
corresponding homomorphism. In other words 

(11.1) tv(x) : T^-^iu) ^ (u - x^R^-^iu - x) . 
Introduce a map V( x i> •• - ,x r ) : Y(gl N _ 1 ) — > End (W® r ) ® ^(fltjv) > 

(11.2) . . . , x r ) = (tt(xi) (8) ... (8) vr(a: r ) <g> V) ° (A <JV_1) ) (r) , 

where (A^ 7V_1 ^) < ' r ' ) : F(flt/v-i) — * (^(f^AT-i)) 8 r+1 i s the multiple coproduct. For any 
XeEnd(W) set v(X) = Iwi . Define a map ^(xi,...,x r ) : ^(fl^-i) -> W® r ® Y(0y, 

V»(xi, . . . , x r ) = (z/® r <g> id) o ip(xi, . . . , x r ) . 
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Lemma 11.1. iP(x u ...,x r ) (^(flU-i)) C W® r © Y + (gl N ) . 

Proof. The claim follows from the facts that Y + (gl N _ 1 ) is a coideal in Y(gl N _ 1 ) , and the 
vector wi is a singular vector with respect to the action of Y(gl N _ 1 ) . □ 



Consider the embedding : Y(gl N _ 2 ) - Y^l^) : 0(T a <f 2> (w)) = T^"^) , 

a, 6 = 1, . . . , AT — 2 , similar to the embedding -0 . Recall that 6 : ^(gt/v) — > ^(gt/vV^+G^Ar) 
is the canonical projection. 

Lemma 11.2. We have (id® r © 0) o ^(xi, . . . , x r ) o = wf r © (0 o ip o 0) . 

For any element # G W 8 ""© Y(gl N ) we define its components ^r a iv,«r ^fie ru ^ e 

JV-l 

9= E w ai ©...© w Sr ®5 ai ar . 

Set e= (e 2 e^- 1 ) and t = (f?, . . . , t£ a ; . . . ; t? ~\ . . . , tp_\) . 

Proposition 11.3. [TV1, Theorem 3.4.2] We have 

N-1 

(11.3) B € (t) = £ T^^)...^,^^ 

ai,...,a ? i = 1 

Proof of Lemma 5.2. Lemma 11.2, Proposition 11.3 and formulae (11.1), (11.2) imply that 
the denominator of 9(M^(t)) is at most Yl Ti Yi tf) • Now the statement follows 

a=l i=l j = l 

from formula (5.2). □ 

Since V = Cvi © W we have V Ak = (vi A W A(fc_1) ) © W Afc , where the first summand is 
spanned by the vectors Vi A v ai A ... A v afe _ 1 , 2 ^ a\ < . . . < ak-i ^ N , and the second 
summand is spanned by the vectors v^, 1 A ... A v^ fc , 2 ^ bi < . . . < bk ^ N . Abusing 
notation, we will often identify W A ^ _1 - ) and vi A W A( ^ fc_1 - ) taking a vector x to vi A x. 

The R- matrix -Raa^aiC^) acting in V Ak ®V preserves the subspaces 

(vi A W A(fc " 1) )©Cvi, ((viA W A(fc - 1} ) © W) © (W Afc ©Cvi), W Afc ©W. 
Its restrictions onto these subspaces are 

(11.4) «Afc,Ai(«)| (viAWrA(fc _ 1))0Cvi = « + J R Afc , A1 («)|^ H/ = 4 i I- 1 1 ) («), 



Rf\k,/\l(u) 



((viAVK A(fe - 1) )(g)W)©(W Afe (g)Cv 1 ) 



K A(fc-l),AlW 15 



where ^((vi Ax)© w) = (w A x) © vi , and we will not use the explicit form of 5" . A 
similar decomposition can be written for the .R-matrix i? A i,Afc(w) . 
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In what follows we regard T(u) and T Ak (u) as matrices over the algebra Y(gl N ) and 
introduce their submatrices induced by the decompositions V = Cvi © W and 
V Ak = (vi A II^" 1 )) © W Afc : 

A{u) B{u)\ mAfc , A fl(u) 



(11.5) Tin) = ' , T Afc (w) = , - 

\C(«) %)/ 1 ' \c{u) D{u) 

For example, A{u) = T n (u) , B(u) = (T 12 (u), . . . ,T liV (u)) , £>(u) = (T lJ (w))^ =2 . 

Denote by HY(L,M) the space Horn (L, M) © Y(gl N ) of matrices with noncommut- 
ing entries. We call L the domain of those matrices. The coefficients of the series 
B,D,A,B,D respectively belong to HY(W,C), HY(W, W) , HY(W A ( fc " 1 ), W rA ( fc " 1 )) , 
HY(W Afe , W^" 1 )) , HY(VF Afc , W Afc ) . For k = 1 we have ll^" 1 ) = C , A = A, B = B, 
D = D . 

Set 

Define a map S : W A( - k ~^ © W -> W Afc by the rule 

(11.6) S(x© w) = wA x. 

We will use the following commutation relations that are obtained from formulae (4.11), 
(11.4) , Corollary 3.3 and the equality S R(u - 1) = (u - k) S : 

(11.7) bW(u)bW(v) = U ~ V bW(v)bW(u)rM(u-v), 

u — V + 1 

(11.8) AM(u)bW(v) = B^(v)A^(u)R (12 \u-v-l) + B (u) A{v) , 

(11.9) DW^B^W = B^(v)D^(u)R^ 12 \u-v) — SB^(u) D^(v) . 



u — V 



The equalities hold respectively in RY(W Ak © W, W A{ - k ~^) , BY(W A ^-^ © W, W A ( fc " 1 )) 
and HY(lF Afc © ll 7 , W Ak ) , and the superscripts in brackets indicate what tensor factors are 
domains of the respective matrices. If k = 1 , relations (11.7) - (11.9) become 



(11.10) bW(u)bW(v) = u _ v + 1 bW{v)bW(u) R( 12 \u-v), 

ii — ii — i i 

(11.11) A(u)B(v) = - — - — - B(v)A(u) + — — B(u)A(v) , 

u — v u — V 

(11.12) D^\u)B^ 2 \v) = B^\v)D^(u)R {12) (u-v) — B^(u) D (2 \v) 



u — V 
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Let R(u) = (u + l) -1 P^ N ~^ R^ N ~^(u) . For an expression f(ui,...,u r ) with coeffi- 
cients in matrices with the domain W® r and a simple transposition (z, , i = 1 . . . k — 1 , 
set 

(11.13) (t ' t+1) /(«l,---,Ur) = /(til,..., 'Ui_l,'Ui + l,Wi,Wi + 2,...,'U r ) J R (j ' J+1) (w i -'Ui + i). 

The matrix .R(it) has the properties R(u)R(—u) = 1 and 

^ 12 )(w-z;) J R( 23 )( W )^ 12 )(z;) = ^ 23 )(t;) J R( 12 )( W )^ 23 )( W -t;), 
cf. (3.2). This yields the following lemma. 

Lemma 11.4. Formula (11.13) extends to the action of the symmetric group S r on expres- 
sions f(ui, . . . , u r ) with coefficients in matrices with the domain W® r : f i— > a f , a G S r . 

By formula (11.10) the expression B^(ui) . . . B^ r \u r ) is invariant under the action (11.13) 
of the symmetric group S r . 

For an expression f(u\, . . . , u r ) with coefficients in matrices with the domain W® r , set 

(n.14) ^Sym^;;^ /(«!, ...,u r ) = J2 7(«i, • • • , «r) • 

Proposition 11.5. 

(11.15) l(°)(«)B[ 1 ](M 1 )... J BM(M r ) = 

= Bl 1 ^!) . . . BM(t( r ) l°'(«)^ 0r '(ti - « r - 1) . . . I ^ - «n - 1) + 

+ BMVtlf— f[ Ul ~ Ui ~ 1 S^B^(u 2 )...B^ {Ur) A( Ul )), 

(r-lj! ' r V w -^i ^i - ^ / 

(11.16) 5W(w)S[ 1 ](w 1 )...sM(w r ) = 

= S^^-ui) . . . B^ r \u r ) D(°\u) R^ 0r \u — u r ) . . . 

x R {lr) ( Ul -u r )... R {12 \ Ul - u 2 )j , 

where the tensor factors are counted by 0, 1 , . . . , r . 

Proof. The statement follows from relations (11.7) - (11.12) by induction with respect to r . 
We apply formula (11.8) or (11.9) to the product of the first factors in the left hand side 
and then use the induction assumption. □ 

Remark. Formulae (11.15) and (11.16) have the following structure. The first term in the 
right hand side comes from repeated usage of the first term in the right hand side of relation 
(11.8) or (11.9), respectively. The second term, involving symmetrization, is effectively 
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determined by the fact that the whole expression in the right hand side is regular at u = Ui 
for any i = 1, . . . , r , and is invariant with respect to action (11.13) of the symmetric group 
S r . The symmetrized expression is obtained by applying once the second term in the right 
hand side of the relevant relation (11.8) or (11.9) followed by repeated usage of the first 
term of the respective relation. 

N . 

Let Q = QaEaa G End (V) be the diagonal matrix from Theorem 5.4. Set 

o=l 

N-l 

(11.17) Q=J2 Qa+i ^ia _1> G End (W) . 

a=l 

Denote Q = Q A(fc_1) and Q = Q Ak . Then formula (4.13) implies that 

(11.18) ?k,Q(u) = Qi ti w «k-i)(QA(u)) + tr W Ak(QD(u)) . 

To prove formula (5.8) we will employ Propositions 11.3 and 11.5. In our present notation, 
relation (11.3) reads as follows: 

(11.19) B € (t) = B^(tl)...B^\tl 1 )^tl...^ e )(Mf- 1) (i)). 

We substitute (11.18) and (11.19) into the left hand side of formula (5.8) and use relations 
(11.15), (11.16) to move A(u) and D(u) through the product B^(t\) . . . B [€ll (*Ji) . 
To simplify writing, let r = and Ui = t\ , % = 1, . . . , r . We have 

(11.20) <J k , Q {u)B^\u 1 )...B^{u r ) = 
= B^( Ul )...B^(u r ) x 

x (Qi(tr WA(fc -u <g> id® r ) (g(°) (u) R^ r \u - u r - 1) . . . # 01 >(u - Ul - 1)) + 

+ (tr WAfc <g) id® r ) (Q (0) 5 (0) (u) £ (0r) (u - u r ) . . . £ (01) (u - ui))) + 
— % l £(u)BW(u2)...BW(u r ) x 

U — U\ v 

fai 17 A( Ul ) - QWDW(tt 1 )F^« 1 -« r )...F 2 )(« 1 -« 2 )) , 

V - M < /J 

where the tensor factors for the products under the traces are counted by 0, 1, . . . , r . Here 

(11.21) S Q (u) = (tr w/A(fe - 1) ®id)(g A ( fc - 1 )%)5) 



x 
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is a series with coefficients in HY(VF, C) . Besides this, we use the equality 

(tr w , k ®id)(Q Ak S(B(u)®idf) = S g (w)Q 

that follows from the cyclic property of the trace and the formula Q Ah S = S (<5 A( - fc_1 - ) <E> Q) ■ 
Notice that Bq(u) = B(u) if k = 1 . 

To proceed further, we use the next lemmas. 

Lemma 11.6. 

It — U r ) - - - 

£(01) 

(it — tti) = ip(ui, . . . , u r )(T( N ^(it)) , 

D(°)( U )fl( 0r )( U - Ur )...^ 01 )(u- Ul ) = V(«i,...,«r)((T< iV - 1 >(«)) Afe ). 

Proof. The equalities follow from formula (4.5), Lemma 4.3 and formula (11.2). □ 
Lemma 11.7. For any X G ^(fltjv-i) > 

(11.22) A(u)ip(u 1 ,...,u r )(X) ~ V(ui,...,u r )(X)A(u), 

(11.23) A(°\u)R( 0r \u-u r -l)...R i01 \u-u 1 -l)ifj(u 1 ,...,u r )(X) ~ 

- n aw . . . , «r) ((T^-^n - i)) A(fc - i} x) . 

Proof. Denote by Y^gl^) the left ideal in Y^gl^) generated by the coefficients of the 
series T a i{u) for 2 ^ a ^ N . It is a subideal of Y + (gl N ) . Relations (4.2) imply that for 
any Z G Y(gl N _ 1 ) and CeY x (gl N ) the coefficients of the series [T-q (it) ,^(Z)] and the 

product CV>(Z) belong to ^(flU). Then formula (11.22) holds because A(u) =T^\u) 
and the coefficients of components of ip(ui, . . . , u r )(X) belong to ip(Y(gl N _ 1 )) . 

To prove formula (11.23) we first observe that the coefficiens of entries of the matrix 
A(u) — A(u) _D A ( fc_1 ) (it — 1) belong to Y x (gl N ) ; this fact follows from the formula for A(u) 
similar to the first formula of (4.12) for the identity permutation a . Then we proceed like 
in Lemma 11.6 and use the fact that ip(ui, . . . , it r ) is a homomorphism of algebras. The 
calculation goes as follows: 

(u — ui — 1) ^(iti, • • • , u r )(X) ~ 
~ A(u) (L> A(fc - 1} (ix - 1)) [ ° ] R {0r \u -u r -l)... R {01 \u - Ul - 1) • • • , u r ){X) ~ 

~ A(u) </>(«!, • • • , U r ) ((T^iu - l))^) • • • , « r )(X) TT ~ 

1=1 

* A{u) if>( Ul , ...,u r ) ((T< N -»(u - l)) A(fe - 1} x) n ■ 

. = 1 U Ui □ 
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For an expression f(v) set res f(u) = ((v — w)f(v))\ , if the substitution makes 

v — w \v = w 

sense. Observe that 

(11.24) Q^D^im) R {lr) ( Ul - u r ) . . . R {l2 \ Ul - u 2 ) = 

= res ((tr w <g> id® r ) (Q (°) £>(°) (v) R {0r) ( i> — w r . . . 



because res (R(v — ui)) =p( N - 1 ), 

V — Ul 

Altogether, for any X e ^(flt/v-i) we have 
(11.25) T fciQ (u) . . . B [r] (u r ) 1>( Ul , u r ){X) ~ 



B^K) . . . BW^) ...,«,.) pr£-g(« - 1) X) Q t A(u) J] 



M — Ui — 1 



+ 



U — Ui 

i=l 1 



+ ^ 1) ..., Ur )(T^-V)X)) + 



1 i? Q (l-r) 



— l B l Q ] (u)B^(u 2 )...B^(u r ) x 



M - Ml Q 



x U(u u . . . , u r )(X) Qi n U \ 1 - res^wi, ...,«,.) (O^^t/) X))) 1 } . 

The next step is to apply both sides of formula (11.25) to the vector wf r . That 
amounts to replacing there the map ip(ui, . . . , u r ) by ip(ui, . . . , u r ) , and changing the sym- 
metrization ^Sym^ r Ur , cf. (11.14) , to the ordinary symmetrization Sym Ui Ur , because 
R(u) wi <g> wi = wi (g) wi . 

Now we replace X in formula (11.25) by B^ -1 ^?) , cf. (11.19) , and employ the induction 

assumption. We recall that A(u) = Tn (u) , r = , and Ui = t\, i = 1, . . . , r . We also use 
Lemma 4.2 and equalities 

it — X 

tt(x) {T^- l \u)) wi = wi , a = 2, . . . , AT — 1 , 
cf. (11.1), which imply that 



^(^,...,4 1 )(r 1 < f" 1> («)) ^ wf®r 22 («) n 



M - ^ + 1 

u-t] ' 



V(f i , • • • , tji) (Tjf -^H) ~ wf r ® T a+M+1 («) , a = 2, . . . , AT — 1 . 

35 



Besides this, we use Lemmas 4.1 and 11.1. As a result, we transform formula (11.25) to 
(11.26) 7 kyQ (u)M^(t) ~ 

^ i=l 1 a j=2 

k . 

+ E II X ^ "J + + U € , fc>g («;t) , 

6 j = l / 

cf. (5.7), where the first sum is taken over all (k — l)-tuples a = (a,2, ■ ■ ■ , at) such that 
2 < a 2 < • • • < afc < AT, the second sum is taken over all /c-tuples b = (pi, . . . , bk) such 
that 2 ^ 6i < . . . < b k ^ N , and 



(H-27) U c>fc>g («;t) ~ 



* 5W(tD . . . bw^) . . . , 40 (uf fc :;> Q (« - i;o) Qi n ^^r 1 + 



1=1 * 



1 



+ 77\ TTf Sym t i . ,i 



^^\u)B^{t\)...B^{t\,) x 



u — 1 1 



x (^(t!,...,fJ,)(Bf" 1> (f)) x 



V *=2 J l ^ 4 =2 J l 7 = 1 J l 
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- res^,..., t\.) (ug-£ (v;t))) 



The last formula determines II ^ k q(u;t) modulo terms with coefficients belonging to 
Y + (gl N ) . By the induction assumption the expressions U^~^\v;i) belong to I^[~^(v;i) . 
Then together with Lemmas 3.2, 4.1, 4.2, 11.1, 11.2, and formulae (5.5), (11.1) , (11.2), the 
induction assumption implies that there exists an expression IL k q(u; t) which belongs to 
I^,k,Q and obeys formula (11.27). Therefore, formula (11.26) coincides with formula (5.8). 
Theorem 5.4 is proved. □ 

Proof of Proposition 6.2. We prove the statement by induction with respect to N using 
Proposition 11.3. It is straightforward to check that the induction assumption yields 

e ab B«i®-®«« (t;z) = , 2 ^ a < b < N. 
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Let C(u) be the left bottom block in the decomposition (11.5) of T(u) . Consider the 

oo 

coefficients of the series C{u) , C(u) = C s u~ s . Then C\ is a matrix from HY(C, W) 



s=l 



whose entries are ei2, • • • , eiN . Therefore, it remains to check that CiB^ 1 n (t; z) = . 
Relations (4.2) implies that 

dB(u) - B(u)d = D{u) - A{u) . 

Then, similarly to the proof of Proposition 11.5, we get 

CiS [1] (ui)...S [r] (u r ) = B [11 K)...B M (« r )C 1 + 

+ (^ fl Sym^(5[ 2 ](. 2 )...5HK) x 

x (dWk)^^-^...^ 2 '^-^- n Ul ~ u ^ 1 A( Ul ))). 



The rest of the proof of Proposition 6.2 is similar to considerations in the proof of Theo- 
rem 5.4 after Lemma 11.6 and in the proof of Theorem 6.1. □ 

12. Proof of Theorem 8.6 

We will prove Theorem 8.6 by induction with respect to N, suitably adapting the proof 
of Theorem 5.4. 

We will continue to use notation introduced in the previous section. The maps ip , ip(xi, 
. . . , x r ) and ip(xi, . . . , x r ) defined there for the Yangian Y(gl N _ 1 ) have their counterparts 
for the current Lie algebra g[ JV _ 1 [x] . Abusing notation we will denote them by the same 
letters. 

The embedding ip : U(gl N _ 1 [x]) — > U(gl N [x]) is defined by the rule ^(L^ -1 ^)) = 

L a+i) 5 + i('u) , a, b = 1, . . . , N — 1 , where L^ -1 '(it) is series (8.1) for fl[j V _ 1 [a;] . The ho- 
momorphism ip(xi, . . . , x r ) : U(gl N _ 1 [x]) — > End (W® r ) <g) U(gl N [x]) is given by 

, AT , r i ^ Z?^ -1 ) (5?) 1 ®(r-i) 

• • • , x P ) (L^- 1} («)) = 1^® L«+i lfc+ i(«) + £ - ® 1 , 

a, 6= l,...,iV- 1. The map rp(x lt . . . ,x r ) : t/foU-ifc]) -»• ® ^(flU^l) is such 

that 

r _ ®(r-i) 

V»(xi,...,x r )(L^ >)) = wf r ®L a+1>6+1 («) + (5 al 2^ x - 1, 

. . U Xi 

1=1 

a, b = 1, . . . , N — 1 , and V>(xi, • • • > ^vX^i-X^) = V^i) • • • > ^r)(^i) ^(^l, • • • , x r )(X 2 ) for 
any X 1 , X 2 G *7(flljv_i[a;]) • 
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Lemma 12.1. We have 

^(xi,...,x r )(?7(flU_ 1 [x])nf- 1> [ a ;]) C W® r ® U^i^x]) nf _1> [x] . 

The proof is straightforward. 

Define the embedding : £7(0^-2^]) ^ ^(flU-ikl) , ^C^f = L J+i~!b+i( u ) > 

a, 6 = 1, . . . , AT — 2 , similar to the embedding ^ . 

Lemma 12.2. We iiave ip(xi, . . . , x r ) o <p = w® r ® (ip o (f>) . 
The proof is straightforward. 

For any element g G W® 7 "® C/(0ljv— 1[^]) we consider its components g a i>—' a r given by 

N-1 

9= E w 0l <g>...<g> w ar ®j ai flr . 

Recall that £ = (£ 2 , • • • , and f = (t?,...,t| 2 ; ... ; tf "\ . . . , if^) . 

Proposition 12.3. We have 

iV-1 

(12.1) F € (t) = £ L^^)...^^^ 

ai,...,a ? i = 1 

Proof. The statement follows from formula (8.7). □ 

In Section 10 we introduced the degree function deg on the algebra Y(gl N ) [[it , d u , £]] 
by the rule deg u~ x = deg d u = deg ( = — 1 . We define the degree to rational expressions in 
t\, . . . , t^~-\ with coefficients in Y(qI n )[[u~ 1 , d u , Q] , assuming in addition that deg tf = 1 

and deg(t" — t^) _1 = — 1 for all possible a,b,i,j. We use similar definition of the degree 
for the Yangian Y(gl N _ 1 ) and related rational expressions in t%, ■ ■ ■ , t^l\ . Notice that the 
maps ip , ip(t\, . . . ,tli) and ip(t\, . . . , th) agree with the introduced degree functions and 
the projections to the graded algebras. For example, if X G Y(qI n _ 1 ) and deg X = k , then 
degip(X) = k and ip (grad k X) = grad fc (V>P0) . 

Proposition 12.4. We have deg(B c (t)) = , and grad_ e (B^(t)) = F c (t) . 

Proof. The statement follows from formula (10.2) and Propositions 11.3, 12.3 by induction 
with respect to N. □ 

Proof of Theorem 8.6. As in Section 10, set (5 = 1 + C-^- We will employ Theorems 10.1, 
10.2, and Proposition 12.4, and adapt arguments used in the proof of Theorem 5.4 to prove 
the statement of Theorem 8.6 by induction with respect to N. 

Consider formula (11.25) and write there H>k,Q(u) instead of 'Bq(u) to indicate the 
dependence on k . Recall that we have r = itj = t\ , i = 1, . . . , and -A (it) = Tn(it) . 

(— lV~ fc (jV — k)\ 

Taking the sum of equalities (11.25) for k = 1, . . . ,1 , with coefficients — T)V{1 — fc)T 
and using formula (10.3), for any X G ^(fltjv-i) > we obtain 
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(12.2) S l>g («) B^(t{) . . . B^\t\.) if,(t\, . . . ,t\,){X) ~ 
~ ^\t\) . . . B^\t\.) x 

pi 

x (^,...,^)(s^-J>(«-i)x) (gxTuH n 

+ , tJO ((8<?> - 1) - C>) + X)) + 



+ 771 777 ^Y 111 * 1 t 1 



« -fj " ' X 



+ (l - raj • • • , f JO (^% _1> («)))) iK* 1 , • • • , *J0 (*)) ] } 



where 

(12-3) = £ (-I)'"* ^-^ M«) • 

Notice that in formula (11.25) for k = N we have T ) [^" 1 '(m) = and Hn,q(u) = 0, and 

in formula (12.2) for I = N we assume that Sj^^(it) = 0. We remind that To,q(w) = 1 

and T^ _1 ^(u) = S^ -1 ^(u) = 1 by convention. 

Recall that ^B/ C; q( , u) is an (N — 1) x 1 matrix whose entries are series in u~ l with coef- 
ficients in Y(gl N ) , and 7^q(u) is an (N — 1) x 1 matrix of the same kind, 

Bfc,g(u) = (Bfc,Q;l(w),-..,Sfc,Q;iV-l(«)) , ?/,q(«) = (3 , I,0;l(w),---,3 > I,Q;iV-l(«)) ■ 

Lemma 12.5. For any m = k, . . . , N , distinct i\, . . . , e {1, . . . , to} , and any s = 1, . . . , 

fe , we nave 

to! {N -m)\ 
® k ' Q ' a{u) = (k-l)l(N-k)! X 

x (tr v8m ® id)(4+i,i Q in) - ■ ■ Q (tk) T^ m+1 \u) . . . T ^ m+1 \u - k + 1) . 

Proof. Formulae (11.21), (11.17), (11.5), and (11.6) give that 

<B k , Q ., a {u) = k (tr^ fc ® id) (-Eo+1,1 • • Q (fc) T^+^tx) . . . T( k > k+1 \u - k + 1) A<iy *>) . 
The remaining consideration is similar to the proof of Lemma 4.8. □ 
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Lemma 12.5 and formula (12.3) allow us to obtain "Bk,Q- a (u) and J^Q^it) as coefficients 
of a formal difference operator similarly to formulae (4.14), (10.6), 

(12.4) ± (-1)<*-D S fc>g;a („) e (-*)«- = 



fe=i 



(m — k) ! 
(iV-/)! 



= E J^Tfy, *i.Q-M ( e °" - 1 ) m ~' = 

= m(N-m)\(tr v ^ m ^id)( K E^ hl QWT^ m + 1 \u)e- d - x 

x (1 -Q^T^ 2 ' m+1 \u)e- du ) ... (1 - Q^T^ m+1 \u) e~ du ) A^j" ) e m9 « . 
Introduce the following (AT — 1) x 1 matrices 

^.JfH = {3 r i,K;l(u),...,3 r i,K;N-l(u)) , I = 1, . . . , N , 



whose entries are defined by the rule 

JV 

\ 



(12-5) ^J2(-l) (l - 1} ^K;a(u)d u N - l = 

1=1 

= (t W ®id)(i£^^ Af-p) . 

It is easy to see that 3 r i,i<:(w) = 0. Like in the proof of Theorem 10.2, formulae (12.4) 
and (12.5) imply that deg (^^(-u)) = 1 — / for any I = 1, . . . , N, and 

grad 1 _;(T z , Q (w)) = ?ixk{u). 

N-l 

Set K= £ K a+1 £ja \ so that Q = 1 + (K • 

a=l 

Lemma 12.6. We have deg (l - res (V>(*i, • • • , tji) ( T 1 %~ 1> (^)))) = - 1 aild 
grad^ (l - • • • , ^) (^q" 1 ^)))) = 

Proof. The statement follows from formulae (11.24) and (8.5) by straightforward calcula- 
tions. Notice that the denominator of the expression (v — t}) _1 — ip(t\, . . . , t^i) ^C 1 ')) 

does not actually contain the factor t> — £;[ , which makes the substitution v = t\ in the 
expression legal. □ 
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The Poincare-Birkhoff-Witt theorem for the Yangian Y(gl N ) , see [MNO, Corollary 1.23] , 
[Mo, Theorem 2.6] , yields that if X e Y + (gl N ) , deg X = d , then grad d X e U(gl N [x]) n + [x] . 

Let F(u) = (Li 2 (w) , . . . , Lin(u)) = grad_ x (B(u)) , and d = degX in formula (12.2). 
We apply the map grad d _ z to both sides of (12.2), and then evaluate the result at £ = 1 . 
Recall that we have 

grad_ ; (§i,q(u)) = Gi, CK (u) , grad_ s (§^ _1> (w)) = G^^iu) . 

Thus we obtain that for any Z e U(gl N _ 1 [x]) , 

(12.6) g hK {u) F^(t\) . . . Fte\t\i) t l e ){z) ^ 

- Wl) ■ ■ ■ L(t{, . . . , tfr (G^(u) Z) [K, + L U („) - £ ^t) + 

+ ^L....40((e, < r , («)-C>)) z )) + 



' 1 —?\%{u)FW(tl)...FKXt\>)y. 



+ !T\ TTT Sym t i t i 



u-t[ 



x (^{t\, . . . ^{Z) (K x + L xx {t\) -J2 + 



The rest of the proof of Theorem 8.6 is similar to the proof of Theorem 5.4 after for- 
mula (11.25). We apply both sides of formula (12.6) to the vector w® r , that amounts 



• "l J ". N < > " ! (I,.; , fUHL I I U ll l rtM z. — J ' ' 

Since formula (12.1) is equivalent to 



to replacing there the map ip(ui, . . . , u r ) by ^(^l) • • •■> u r) > and then take Z = F| (t) 



F*(f) = F^ft) . . . F^l(tJO • • • , * JO (Ff "^(t)) , 

the left hand side of (12.6) becomes Gi,k{u) F^(£) , and to transform the right hand side to 
the required form, we employ the induction assumption. □ 
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Appendix A. The Gaudin model as a limit of the XXX-type model 

Let Mi, . . . , M n be jjt^-modules. In Sections 6 and 7, we describe the XXX-type model 
on the tensor product Mi ® . . . <g) M n considering the action of the Yangian Y(gl N ) on the 
tensor product Mi(zi) (g) . . . <g) M n (z n ) of evaluation modules. The objects associated with 
the XXX-type model are the operators 

r„r- 8M "M = ^M| MiM8 8MnM , „,» = i iv, 

rationally depending on w and z\ , . . . , z n ; the transfer matrices J^-® M ™ ( u; z ) , c f. (6.1) ; 
the modified transfer matrices 

fc '3 V ' ^ fc ' VV ; lM 1 ( 2l )®...®M„( 2 „) 



1 \^ / -iNfc-i (X- /)! T M 1 ®...®M„/ n 

(iV _fc)! ^ l iJ {k-l)\ J ^ W 



/c = 0, . . . , X, cf. (10.3) ; the difference operator 

N 

(A.l) ®^®-® M "(u;d u ;z) = ]T (-l) k "® M " («; z) e"* - = 

fc=0 

= ^ (_i)* sJ£®-® M » ( u; z) (e a " - l)(^- fc ) 5 « e"^ ; 



fc=0 



the deformed Shapovalov form <Sj^ ,g M on Mi <S> ■ ■ ■ <S> M n , see (7.9); the universal 
weight function W^®'"® Vn {t; z) , cf. (6.2); the rational functions 

1 1 +a ~ . i A a +! 11 +a j.o-1 11 w_wj_l 11 



Va+1 j=1 t, Z 3 + Aj- j=1 t l t J j = 1 l t Zj + 1 j=1 t t tj 1 

a = 1, . . . , X — 1 , i = 1, . . . , £ a , that define the Bethe ansatz equations (6.3) ; and the 
difference operator 

m^ Q (u,d u ;t;z;A) = (l - Xj iQ (w; t; z; A) e" a ») . . . (l - X% Q (u; t; z; A) e~ d ") , 

cf. (6.6) , (6.4) , whose coefficients are eigenvalues of the transfer matrices. 

In Section 9, we describe the Gaudin model on the tensor product Mi <g> . . . <g> M n consid- 
ering the action of the current algebra fll/vM on the tensor product M\(zi) <S> ■ ■ ■ <S> M n (z n ) 
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of evaluation modules over flljvt^] • For the Gaudin model, the associated objects are the 
operators 

n (i) 

L M^M n{u . z) = a=l,...,N, 



rationally depending on u and z\, . . . , z n ; the transfer matrices £^®---® M ™ (u; z) , cf. (9.1) . 
(9.2) ; the differential operator 



N 

(A.2) V^-^ M -(u,d u ;z) = ^2(-l) k g?£"'* Mn (u;z)d?- k 

k=0 



with End (Mi <g> . . . <g> M n )-valued rational coefficients, cf. (8.3); the Shapovalov form S^i® 
S Mn on Mi <g> . . . <g> M n ; the universal weight function F^ 1 ®"-^" (£; z) , cf. (9.6); the 
rational functions 



%lf K (t;z;A) = 



a«-a" +i * i * i t 



K a K a+1 + ^2 3 fa_t. + E +a ,a-l 2 E /a _ /a + E +a ,a+l> 
j = l ' j=l ~ './ j=l ' I j=l l i ~ '.i 

a = 1, . . . , N — 1 , i = 1, . . . , £ a , that define the Bethe ansatz equations (9.3) ; and the 
differential operator 

£CL— 1 AO 

— > / " A« ^ 1 « 1 \ 

1^JV V i=l " i=l M ~ ^ *=1 U l *' 

cf. (9.4) , whose coefficients are eigenvalues of the transfer matrices. 

The objects associated with the Gaudin model can be obtained from the corresponding 
objects for the XXX-type model in the following limit. 

Set e~ x z = (e~ 1 zi, . . . , e~ 1 z n ) and e^t = (e -1 ^, . . . , e~ 1 t^ N ~_\) . 
Theorem A.l. Let Q = 1 + eK . As e — > , we have 

(A.3) 2^®-"® M »(e- 1 «;e- 1 z) = S ab + sL^®-® M "(u; z) + 0(e 2 ) , 

a,b=l,...,N, 

(A.4) S^®"-®^- 1 ^- 1 *)) = e k Q*f£-® M "{u;z) + 0(e k+1 ) , 

k = 0, ...,N, 

(A.5) D^Q'" <s>Mn (e~ 1 u, ed u ; e~ x z) = e N T>^®-® M "(u; d u ; z) + 0{e N+1 ) , 
(A.6) s£l..®M n ~* S Ml ®...®S Mn , 
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(A.7) 



(e-H;e~h) II II II 



n JV-1 C N-2 C C + 1 

£ -n- -n- -n- £ 



nnn^ 



_ 7 . J. J. J. J. J. J. j. a -l-i _ +a 
i=l a =l j = l 3 * % o=l i=l j=l l 3 l i 



(A.8) Q^Q(£~ 1 t;£~ 1 z; A) = 1 + sX%f K (t; z; A) + 0(e 2 ) , 



(A.9) m^Qie-^ed^e-H-^-^A) = e N M^ K (u, d u ; t; z; A) . 

Proof. Formulae (A. 3) and (A.8) are straightforward. The proof of formulae (A. 4) and (A. 5) 
is similar to the proof of Theorem 10.2. Formula (A. 6) follows from formula (7.4). For- 
mula (A.7) is similar to the claim of Proposition 12.4. □ 

Remark. Using Theorem A.l one can see that Theorem 7.6 implies Theorem 9.1, and Theo- 
rem 6.1 implies Theorem 9.2 for an isolated solution t = (£[,..., t^lX) of the Bethe ansatz 
equations (9.3). 

Appendix B. Dynamical Hamiltonians as parts of transfer matrices 

N 

Let d= J2 e aa , and C 2 = ( e a a ebb - e ab e ba + e aa ) . They are central elements 

a=l l^a<b^N 

of the universal enveloping algebra U(gl N ) . 

N 

Let Mi,...,M n be gl^-modules. Let K be a diagonal matrix, K = K a E aa . The 

a=l 

Gaudin Hamiltonians Hi^k(z), ■ ■ ■ , H U: k{z) , see [G], the rational dynamical Hamiltonians 
Gi,k(z), ■ ■ ■ , Gn,k(z) i see [FMTV], [TV2], and the trigonometric dynamical Hamiltoni- 
ans X^k(z), • • • , X NjK (z) , see [TV2] , [TV3] , acting in Mi <g> . . . <g> M n , are the following 
operators: 

N N n Ji)Jj) 

H iiK ( Z ) = E^ e - + E E 



e ab e ba 



a=l a,b=l j=l Z% Zj 

3=£i 



N 

€-ab €-ba &aa 



C*M = E*42 + E Ka _ Kb 

i=l 6=1 a 

b^a 



( \ 2 n K 

^w = -^- + E^ e 2 + E E ^^ + E^V(^ e ^- e -)- 



i=l 6=1 l^i<j^n 6=1 

b^a 

n . 
W 



Recall that e ab acts on M\ ® . . . <g> M n as e ab 

i=i 
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It is well known that the Gaudin Hamiltonians come from the residues of the transfer 
matrix #-® M »(u; z) : 



gM 1 ®...®M r< 



'2,K 



(u;z) = 



- ^ + E ^ tf> « + E ^ - + E 

1=1 L x 7 = 1 J 7 J 1 = 1 



The rational and trigonometric dynamical Hamiltonians also can be recovered from the 
transfer matrices 0j^?-® M "(u; z) and TjJ^®-® M "(u; z) , respectively. To this end, con- 
sider the following End (Mi <g> . . . <g> M n ) -valued rational function of u, x, z±, . . . , z n : 



G K (u,x;z) = f2(^) k ^f-^(u;z)x N - k f[ -1 



k=0 
N 



a=l 
N 



x - K n 



fc=0 



a=l 



X - K n 



Proposition B.l. We have 



N 



g K (u,x;z) = 1 - u 1 ~Z 



and 



— u 



N 



E^fs-M-ETT^W" 3 * 

o=1 x j\ a v b=i j\ a j\ b j 



K n 



7 K (u,x;z) = 1 - u X J2 



a=l 



— u 



N i<r / ( \ 2 N \ 

E A + ^ - E ) + 0( "~ 3) • 



6=1 

6/o 



Proof. To get the first formula we compute the expansion of G% I }f > "'® Mn (u; z) as u — > oo 
using formulae (8.3), (8.4), (9.1), (9.2), and Proposition 8.1. For the second formula we 
compute the expansion of T^^'"® M " (m; 2) using formulae (4.13), (6.1), and the definition 
of the module M 1 (zi) ® . . . <g> M n (z n ) over the Yangian Y^jj^) . □ 
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Appendix C. Eigenvalues of transfer matrices for real zi,...,z n 

Denote by the real part of the Lie algebra gi N generated by e a b , a,b = 1, . . . , N . 
Let Mi,...,M n be finite-dimensional irreducible gl^-modules with highest weights Ai, 
. . . , A n , and highest weight vectors v\,...,v n . Denote by Mf, . . . , M% be their real parts 
generated over qI^ by the vectors v\, . . . ,v n , respectively. It is clear that for real u, z\, 
. . . , z n , real matrices Q,K, and any k = l,...,N, the operators T^q®"'® m " (m; 2;) and 
gM^®...®M n ^ preserve the tensor product Mf ® R . . . <g) R M* . 

Since Mi,...,M n are finite-dimensional irreducible gl^-modules, the tensor product 
Sm x ® ■ ■ • <8> 5'm„ of Shapovalov forms restricted to M^® R . . . ® R M* is a positive definite 
bilinear form. Hence, an operator, which is symmetric with respect to this form, has real 
eigenvalues. Therefore, Theorems 9.1 and 9.2 yield the following statement. 

Theorem C.l. Let K be a reai diagonal matrix, z\,...,z n real numbers and t = (t\, 
t^ N ~\) a solution of the Bethe ansatz equations (9.3) . If F^ 1 ®"'®'"" (t; z) ^ 0, then for 

any k = 1, . . . , N, the function Z^i-,k(u; t;z;A), cf. (9.5) , takes real values for real values 
of u . 

For a dominant weight A = (A 1 , . . . , A N ) , set A' = A^+ 1 - max{a | A a > A^} . 

For finite-dimensional irreducible gl^-modules L,M, the .R-matrix R LM (u) is a rational 
function of u . Its poles and degeneracy points are known from the representation theory of 
the Yangian Y(qI n ) , see [Mo]. In particular, the following proposition holds. 

Proposition C.2. The R-matrix Rmm( u ) js well-defined and nondegenerate for real 
values of u such that u > A} — A'j or u < A^ — Aj . Moreover, R M M ,(u) has a simple 
pole at u = A[ — A] with the degenerate residue, and the value R M . M .{A\ — A^) is well- 
defined, non-zero, and degenerate. 

It is clear that for real zi,...,z n , the operator R Ml M (z) , cf. (7.6), preserves the 
tensor product Mf ® R . . . <g) R M^ . 

Proposition C.3. Let z\,...,z n be real numbers such that Zi — Zj >A\ — A'- for any 
i j j h3 = l,---,n- Then the bilinear form ®M n i (^-9)> restricted to Mf ® R 

. . . ® R M* is positive definite. 

Proof. The statement follows from Proposition C.2 and formula (7.4). □ 

Theorem C.4. Let Q be a real diagonal matrix, z\, . . . , z n real numbers such that 

Zi — Zj > A] — A'- for any i 7^ j , i,j = 1, . . . , n , and t = (t\, . . . , t^~\ ) a solution of 

the Bethe ansatz equations (6.3) . If (t; z) ^ , then for any k = 1, . . . , N, the 

function 

k 

E l[X^ Q (u-r + l;t;z;A), 

l^ai<...<a fe ^7V r=l 

cf. (6.4) , takes real values for real values of u . 

Proof. The statement follows from Theorems 7.6 and 6.1, and Proposition C.3. □ 
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